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SUMMARY 


A high-performance constant-temperature hot-wire anemometer has been designed based 
on a system theory analysis that can be extended to arbitrary order. A motivating factor 
behind the design was to achieve the highest possible frequency response while ensuring 
overall system stability. Based on these considerations, the design of the circuit and the 
selection of components is discussed in depth. Basic operating instructions are included 
in an operator's guide. The analysis is used to identify operating modes, observed in all 
anemometers, that are misleading, in the sense that the operator can be deceived by in- 
terpreting an erroneous frequency response. Unlike other anemometers, this instrument 
provides front panel access to all the circuit parameters which affect system stability and 
frequency response. Instructions are given on how to identify and avoid these rather sub- 
tle and undesirable operating modes by appropriate adjustment of the controls. Details, 
such as fabrication drawings and a parts list, are provided to enable the instrument to be 
constructed by others. 
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1. BACKGROUND 


A hot-wire anemometer is an electronic device which passes an electric current through 
a fine filament which is exposed to the flow. The filament consists of a material which 
possesses a temperature coefficient of resistivity i.e. as the temperature changes, so also 
does the resistance and hence the Joule heating. The variation of resistance is used to 
generate signals which are related to the flow velocity or temperature. The instrument 
is therefore a thermal transducer which is capable of measuring instantaneous velocities 
and/or temperatures. 

There are two modes of operation of a hot-wire system. The first mode is called the 
constant current mode in which the wire current is kept very nearly constant. The varia- 
tions in wire temperature (and hence wire resistance) caused by the flow are measured by 
monitoring the voltage drop across the filament. The second mode is called the constant 
temperature mode. The filament is placed in a Wheatstone bridge and a feedback amplifier 
is used to maintain the wire at very nearly constant resistance (and hence temperature). 
Fluctuations in the cooling of the filament cause fluctuations in the wire current which lead 
to voltage fluctuations which are measured at the top of the bridge. As the title suggests, 
this report is solely concerned with the constant temperature mode of operation. 

Many workers have the preconceived idea that the hot-wire technique is fully developed 
and that useful measurements can be obtained by simply purchasing a commercially avail- 
able system and following the instruction manuals. After all the instrument has been used 
for more than 50 years and the operating principle appears to be quite simple. However this 
attitude is rather naive and it can often lead to inconsistent measurements or results that 
differ from those obtained by others elsewhere in supposedly identical flows. One source 
for this disparity is that hot-wire measurements are indirect in the sense that variation 
in fluid flow properties are translated into electrical signals via a rather complicated heat 
transfer process which is not yet completely understood. Furthermore, the filament has 
thermal inertia with an associated time constant. The electronic circuitry also has its own 
set of time constants and complicated interactions can occur between the filament and the 
circuitry. In summary, the hot-wire anemometer is a complex instrument and departures 
from correct operation can be insidious and subtle in nature. 

In the late 1960s low-cost high-performance integrated circuit amplifiers became widely 
available and it became relatively simple and inexpensive to construct in-house constant- 
temperature hot-wire anemometers. One example of these 'new’ anemometer designs was 
implemented at the University of Melbourne, based on the ideas of Perry and Morrison 
(1971). The motivation for building these instruments was because of the cost of, and 
the anomalous data obtained with, commercially available instruments. The results of 
various tests and other observations provided convincing evidence that the behavior of 
this particular design was thoroughly understood. In the first year of his Ph.D. (1976) the 
author undertook the task of updating this design e.g. by utilizing up-to-date amplifiers 
which offered higher performance. During the course of this exercise many phenomena 
were observed that were inexplicable in terms of existing theories. These phenomena were 
identified as being of higher than the model of Perry and Morrison. Fortunately they could 
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be recognized during a square- wave test. The operating point could then be changed by 
adjusting the anemometer controls so that these undesirable phenomena could be avoided. 

While at Princeton University during 1986-87 the author worked on methods for in- 
creasing the frequency response of this type of anemometer design for use in supersonic 
flow. A frequency response requirement of 500 kHz is not uncommon in these high speed 
flows e.g. to maintain a spatial resolution of 1mm when the mean velocity is 500 m/s. 
Yet with extremely careful tuning of the system controls by the most highlv-skilled opera- 
tors and under the most favorable of circumstances a barely adequate frequency response 
of around 250kHz can be obtained with commercially available equipment. However the 
so-called phenomena of strain-gauging (in the form of small amplitude high frequency os- 
cillations) was responsible for contamination of about 3 out of every 4 experimental runs 
(A.J. Smits, private communication). Sometimes the oscillations could be made to disap- 
pear by adjusting the anemometer controls. This observation led the author to propose 
that strain-gauging is a purely electronic, rather than an electro-mechanical phenomenon 
and that a more sophisticated model of the system behavior could lead to further under- 
standing and possible control of this frustrating problem as well as leading to methods 
for increasing the frequency response. In order to explore this and other ideas the theory 
of operation of the constant-temperature hot-wire anemometer was extended from 3 rd - 
order to 7 th -order by Watmuff (1987) by including the effects of bridge capacitance and 
the frequency response characteristics of the feedback amplifier. A computer model was 
created based on this analysis. The results of an extensive set of simulations using a wide 
range of parameters led to the conclusion that incorporation of additional controls into the 
anemometer could allow the system frequency response to be increased substantially while 
maintaining system stability. These ideas were not translated into hardware at Princeton. 

In late 1987 the author was sponsored by the Center for Turbulence Research to work 
with R.V. Westphal in the Fluid Mechanics Laboratory at NASA Ames Research Center 
on a project concerning turbulent boundary layers. Westphal had experienced frustra- 
tion with system instabilities when using subminiature hot-wires (e.g. 0.6pm diameter). 
The computer simulation was used to demonstrate that inadequate frequency response 
of the feedback amplifier was the most likely cause for these instabilities. On the basis 
of all of these findings the development of a new high-performance constant-temperature 
anemometer was designed and built which is the subject of this report. 

More recently, a completely general analysis has been performed in which the derivation 
of the system transfer functions have been “algorithmized”. Transfer functions can be 
derived for systems of arbitrary complexity. This work is published here for the first time. 

An effort has been made to include operating instructions. The hot-wire anemometer is 
a difficult instrument to operate properly and detection of an incorrect mode of operation 
can be quite subtle. Examples are given of incorrect operation modes which draw on the 
analysis mentioned above for an explanation. Finally, the report contains a description 
and a schematic of the circuit together with a parts list with suppliers. This information 
is sufficient to enable the anemometer to be constructed by others. 
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2. CIRCUIT BOARD DESIGN 


2.1 Introduction 

Prototype layout arid testing of subsets of the circuit design was carried out by Mr James 
Fogarty and Mr.* Stanley Mason working for Raman Aeronautics. Raman Aeronautics was 
an electronics design and fabrication company contracted to work for NASA Ames Research 
Center over a period of many years. Although Raman Aeronautics is no longer in business. 
Mr. Fogarty and Mr. Mason are presently working at NASA Ames Research Center as 
contractors with Calspan. Testing was also performed by many of the FML personnel, 
in addition to the author. In particular, Mr. Gregory Zilliac performed extensive tests 
on the output amplifier section and was responsible for the discovery of a small 2 nd -order 
nonlinearity. A simplified schematic of the printed circuit board (PCB) and controls is 
shown in figure 2.1. The complete schematic is shown in figure 2.2. 

2.2 Wheatstone Bridge 

The resistors forming the upper half of the Wheatstone Bridge are mounted on the 
PCB and are shown in the upper left hand corner of the detailed schematic in figure 2.2. 
An important design consideration was to minimize the change in bridge resistance with 
temperature as part of an effort to reduce the overall drift of the instrument to a negligible 
level. A relatively large power rating of one watt for the resistors is specified to minimize 
the temperature rise above ambient owing to Joule heating. The resistors are metal film 
with a temperature coefficient of 50ppm per °C so that the drift is even further reduced. 

The standard cross-bridge ratio of 10:1 (i.e. R c /R a = 10) is implemented by inserting 
jumper J 3 while ensuring that jumper J 4 is omitted. The ten lfrU resistors labeled 
RB1...RB10 act in parallel to form the bridge resistor R a — 10011. The resistor labeled 
RB 11 forms the bridge resistor R c = lArfl. Using ten resistors for R a and a single resistor 
of the same type for R c means that the temperature rise of each resistor should be almost 
identical as well as being quite small which will further aid in the reduction of drift . The 
optional symmetric bridge is configured by removing the resistor labeled RB 11, removing 
the jumper J 3 and installing jumper J 4. These operations lead to R a = R c = 20011. In 
the standard 10:1 cross-bridge ratio configuration the current flow through the passive arm 
of the bridge is around 10% of the wire current. With the optional cross-bridge ratio of 
unity the total current capacity of the power supply may need to be increased since the 
current in the passive arm of the bridge will be approximately equal to the wire current. 

The components Rt> and Lt which form the balance arm for the standard 10:1 bridge are 
located on the front panel. A group of resistor decade (1-2-2- 2 code) thumbwheel switches 
is used to set the hot resistance of the probe. The switch is labeled on the sketch of the 
front panel controls (see figure 3.1) and it is configured such that the setting corresponds 
to the balance resistor setting in ohms i.e. the setting shown in figure 3.1 corresponds to 
Rb = 123.011. For the standard cross-bridge ratio of 10:1 this setting corresponds to a 
hot operating resistance of R w ~ 12.31 1 for the wire. The tunable Cable Compensation 
Inductor Lb was selected for use with a co-axial probe cable of total length approximately 
5 meters. The adjustment range is satisfactory for systems with a frequency response 
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FIGURE 2.1 Simplified schematic of circuit. 
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of 20kHz to 30kHz, but smaller values of the inductor may be required for optimum 
frequency response when smaller wires are used since they have a higher frequency response. 
(Refer to section 6.5 for suggestions and examples for improving the stability when using 
subminiature wires). One simple solution in this situation is to use a longer cable rather 
than substitute a smaller inductor. However the overall performance of the system will 
be improved if a smaller inductor is used. In this case the benefits of reducing the cable 
length to a practical minimum should be considered in conjunction with the substitution 
of the inductor since performance is improved when the total inductance is reduced. The 
scaling for selection of a suitable replacement inductor can be based on the ratio of the 
cable lengths. 

The BNC connector labeled EXT RB is used for the optional configuration in which the 
cross-bridge ratio is set to unity. For this optional symmetrical bridge configuration the 
Hot-Wire Resistance Thumbwheel switch and balance inductor must be removed from the 
bridge. The balance resistor Rb should be located at the end of a length of co-axial cable 
that has been carefully chosen to match the probe cable. This configuration offers superior 
compensation to the simple lumped inductor (by using an actual length of cable instead). 
Once again, it is emphasized that the bridge must be symmetrical for this option to work 
properly. 

2.3 Feedback Amplifier 

As mentioned previously many of the ideas for this design were based on the findings of 
the analysis described in section 5. In particular, front panel control would be required for 
the gain and frequency response of the feedback amplifier in addition to the offset voltage. 
In some situations it is possible for the overall frequency response of the anemometer to 
exceed the open-loop frequency response of the feedback amplifier alone. However this is 
not always the case and the feedback amplifier should have the potential for both a high 
gain and a high frequency response in order for the anemometer as a whole to achieve a high 
frequency response. A realistic maximum frequency response for a constant- temperature 
hot-wire anemometer was taken to be 500kHz as could be required in a supersonic flow 
for example. The specification for the feedback amplifier was set at twice that frequency 
(i.e. 1MHz) at the representative gain of A =1000. To achieve such a large gain-bandwidth 
product (i.e. GBW = 10 9 ) it was recognized that the feedback amplifier would have to 
consist of a cascade of individual amplifiers. 

At the outset, the feedback amplifier was designed to consist entirely of operational 
amplifiers conforming to the industry standard 8-pin configuration e.g. LM741. This 
would provide the largest possible choice of devices at present as well as allowing for 
upgrades in the future as higher performance devices become available. The high input 
impedance of instrumentation amplifiers makes them attractive for the first stage of the 
feedback amplifier cascade and a prototype was constructed using three distinct operational 
amplifiers. Unfortunately the noise level was less than satisfactory. The magnitude of the 
noise was found to be dependent on the individual amplifiers even though a variety of 
high quality devices were tried e.g. Analog Devices AD744. Further, the noise level would 
change by up to a factor of five upon substitution of a supposedly identical amplifier i.e. 
one with the same part number from the same manufacturer. Therefore a number of 
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carefully matched integrated instrumentation amplifiers were tested in an effort to reduce 
the noise level. The best combination of performance and low noise was obtained with the 
high performance Burr Brown IN A 103 Instrumentation Amplifier. The advertised gain- 
bandwidth product is 100MHz at a gain K —1000 and the noise level is very low at 1.5 
n\y n/Hz at 1kHz. The instrumentation amplifier is labeled U 1 in figure 2.2. Jumpers J6 
and J1 are used to connect the bridge to the amplifier and they are easily removed for 
testing purposes. The default configuration for {71 is for a gain of 100 which is imposed 
by inserting jumper J5 on the PCB. The amplifier will have unity gain if the jumper is 
omitted. The gain may be adjusted to a nonstandard value by inserting a resistor Rl 
which is selected according to Rl = 6/(K — 1 )k£l. Integrated Circuit numbers U 2, U 3 
and U 4 were originally used for the discrete amplifier components of the instrumentation 
amplifier. These numbers are no longer used since the amplifiers were replaced by the 
single IC Ul. 

The rest of the feedback amplifier cascade consists of four operational amplifiers that 
conform to the industry standard 8-pin layout. The first two amplifiers are AD744 ({75 and 
U 7) and they are configured as variable gain blocks while the other two {U 9 and i712) are 
JFET input LF356 which form the low-pass filter section. The roll-off frequency of each 
variable gain block operational amplifier is adjusted by switching the feedback components 
i.e. the resistors and capacitors which connect the output to the summing junction. The 
roll-off frequency of each filter block is adjusted by switching the resistor which is part 
of an RC network forming the input impedance of the amplifier. One way of performing 
these switching operations is to place the components directly on a front panel switch. 
However this is bad design practice since either long routing distances or long leads would 
be required for the connection from the amplifier to the switch and back. The feedback 
components should be positioned as close as possible to the amplifier for minimum noise 
and for best performance. This objective was achieved in the current design by connecting 
the feedback components with LF1350S analog switches ({76, {78, {710 and U 11) which are 
mounted on the PCB close to each operational amplifier. The analog switches are in turn 
controlled by the front panel switches which generate the digital logic control signals. 

The BCD rotary switch labeled GAIN is used to control the overall gain of the feedback 
amplifier. The BCD outputs drive a logic circuit consisting of the integrated circuits 
Z724, {725 and ?726. The logic circuit controls the analog switches in a manner such that 
the gain of each stage follows the sequence listed in the table on the schematic in figure 
2.2. This arrangement tends to distribute the total gain requirement evenly among the 
two amplifiers. The advertised gain-bandwidth product of the Analog Devices AD744 
operational amplifier is 13MHz so that the frequency response would only just fall short 
of the design specification of 1MHz at the maximum gain level of 16. This corresponds 
to an overall gain K =25,600 for the feedback amplifier. Tunable capacitors are provided 
for the input impedance (e.g. C2) and in the feedback loop (e.g. C6) so that adjustments 
can be made to ensure that the Bode diagram is flat with frequency. The Bode diagram 
(frequency response) cannot be guaranteed to be flat out to these -3dB points unless these 
capacitors have been tuned manually. (See sections 3.4 and 3.5 for tuning instructions). 
The BCD rotary swdtch labeled FILTER is used to control the roll-off frequency of the 
feedback amplifier. Only the first three BCD outputs are used which means that only 
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the first S positions are sequential i.e. switch positions 9 and 10 have the same effect as 
positions 1 and 2. The switch drives the logic inputs of both analog switches directly. The 
approximate -3dB roll-off frequency of the amplifiers follows the logarithmic sequence list ed 
in the table on the schematic in figure 2.2. Note that the LF356 operational amplifiers £79 
and £712 are essentially acting as buffer amplifiers. This arrangement isolates the second 
AD744 gain block from the capacitance of the second filter block. Without this feature the 
capacitive load of the second filter block could cause the AD 744 to undergo high frequency 
oscillations. A simple pole of fixed frequency is incorporated into the feedback loop of each 
LF356 amplifier to ensure stability. The roll-off frequency of the pole is around 2.1 MHz 
which is well beyond the design specification so that it does not interfere with the overall 
performance objective. 

The d.c. offset voltage E qt and the square-wave perturbation voltage are injected 
into the first gain stage. The d.c. offset is adjusted using the front panel potentiometer 
labeled RVA. The square- wave input is derived using the logic circuit consisting of inte- 
grated circuits labeled £730, £731 and £732. A toggle switch is used to provide power to 
the logic circuit and a LED is used to indicate the On position. The circuit is powered 
down when the square- wave is not being used because this helps minimize the noise. The 
amplitude of the square-wave is adjusted using the potentiometer labeled RV 2 which is 
located on the printed circuit board. The square- wave amplitude should be small enough 
to avoid nonlinear behavior i.e. the response to positive and negative fluctuations should 
be symmetrical. A square- wave amplitude of approximately ±100mV has been found to be 
satisfactory. The frequency of the square-wave is adjusted using the rotary switch labeled 
FREQUENCY on the front panel. The square- wave frequency is not a critical parameter 
and control is provided for convenience of observation on an oscilloscope. 

The last LF356 operational amplifier in the feedback amplifier cascade (£712) cannot 
provide the current required to drive the Wheatstone Bridge. The power transistors T1 
and T2 are included for this purpose. The transistors are configured as a Darlington pair 
and testing has indicated that they are capable of supplying a current of 100 mA while 
maintaining a frequency response in excess of 1MHz. The output of £712 is connected to 
the power transistors via the LF13508 analog switch £714 which is toggled by the front 
panel Standby/Operate switch. The arrangement of these components can be seen on the 
second page of the schematic shown in figure 2.2. A LED is used to indicate when the 
switch is in the Operate position in which case the output of £712 is connected to the base 
of Tl. When the switch is in the Standby the position the input to the power transistors is 
grounded and the anemometer is disabled. The switch should be in the Standby position 
when the probe is connected or removed from the instrument otherwise transient currents 
may cause the probe to burn out. The output of the power transistor T 2 is available at 
the front panel BNC connector labeled TOB (Top Of Bridge). 

2.4 Output Amplifier 

Typically, the output voltage of a constant-temperature hot-wire anemometer consists 
of a relatively large d.c. component upon which are superimposed the small fluctuations 
which are of interest. It is common practice to utilize a separate "buck and gain” amplifier 
to subtract a d.c. voltage from the anemometer output and to amplify the resultant signal 
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so that it fills the voltage window of an analog- to-digital converter. For convenience, 
an output amplifier has been incorporated into the instrument for this purpose. The 
output amplifier is essentially an analog computer circuit which is hardwired to perform 
the following linear transformation. 

E 0 = -G(Ei + £ 0 ff) (2.1) 

where G is the gain, E t is an input voltage and E 0 f[ is the offset voltage. It is configured 
as an inverting amplifier because the sensitivity of the anemometer (TOB) output voltage 
is negative. With this arrangement the overall sensitivity of E 0 is positive with respect to 
velocity. 

The front panel Standby/ Operate switch is also used to toggle the LF13508 analog 
switch labeled U13 which connects one of two signals to the output amplifier. When the 
switch is in the Operate position, the input voltage Ej consists of the output of T2 i.e. the 
TOB output voltage of the anemometer. When the switch is in the Standby position the 
input voltage E, consists of the signal applied to the Calibration Input BNC connector 
located on the front panel. The gains of the output amplifier listed in the table in figure 
2.2 are only approximate. Careful matching of the resistors would be required for the gain 
values listed in the table to be realized with precision. However in most situations this 
is unnecessary. The ability to apply known voltages to the Calibration Input allows the 
constants G and E 0 g in equation (2.1) to be determined by direct calibration. One reason 
for wanting these constants with high precision is to enable the recalculation in software of 
the TOB voltages from the digitized values of the output amplifier signals. This would be 
the case for a hot-wire calibration scheme based on King’s Law, for example. However in 
other circumstances it is not even necessary to know the magnitude of the constants since 
the calibration of the output amplifier can be absorbed into the hot-wire calibration as a 
whole. This would be the case with polynomial based calibration schemes, for ex am ple 

The operational amplifier 1734 serves as a voltage follower buffer between the output 
of the LF13508 analog switch 7 13 and the unity-gain inverting amplifier 1715. In the 
prototype design the analog switch was directly connected to the inverting amplifier 1715. 
The input impedance of amplifier 1715 is only 5 k£l and it was discovered that the small 
nonlinearity associated with the impedance of the analog switch introduced a small 2 nd - 
order nonlinearity into equation (2.1). Using a larger input impedance of lOOJbQ for 1715 
reduced the nonlinearity but unfortunately had the side effect of increasing the noise. The 
best solution consisted of the introduction of the buffer U 34 since the extremely high input 
impedance reduced the nonlinearity to a negligible value. This also allowed the low 5kQ 
input impedance to be retained with U15 thereby maintaining the lower noise level. 

The offset voltage E 0 ff is adjusted using the locking counting dial mounted on the panel 
to set the potentiometer labeled RV 3. Unlike the feedback amplifier, the offset and input 
voltages are introduced into the separate unity gain inverting amplifier U 15 rather than 
the first of the variable gain stages. The extra inverting amplifier is available because of the 
requirement of producing an overall negative sensitivity as per equation (2.1). The gain G 
is adjusted using the BCD rotary switch labeled GAIN in the output amplifier section on 
the front panel. The BCD outputs drive a logic circuit consisting of the integrated circuits 
U 27, U 28 and V 29. The logic circuit, the analog switches Ull and 7719 and the amplifiers 
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[’16 and [ T 1S are essentially a duplicate of variable gain blocks of the feedback amplifier. 
The overall gain corresponding to the switch positions follows the logarithmic .sequence 
listed in the table in figure 2.2. Also in a duplicate manner to the feedback amplifier, the 
front panel rotary BCD switch labeled FILTER, the analog switches [’21 and U 23 and the 
operational amplifiers [’20 and [ T 22 serve as simple pole filters for the output amplifier. 
The approximate roll-off frequency of the amplifiers, as defined by the -3dB point, follows 
the logarithmic sequence listed in the table on the schematic in figure 2.2. In an earlier 
version of the circuit the output consisted of amplifier U 22 which was found to undergo 
oscillations when driving long cables. The final stage of the output amplifier now consists 
of an OP-4633 buffer amplifier [’33 which avoids this problem. It should be noted that 
this device is capable of supplying a large current. The lOOfl current limiting resistor R97 
is placed between the output and the BNC connector for protection purposes. In many 
circumstances it may be possible to use the last filter stage amplifier [ 22 for the output 
without any problems and U 33 may be replaced with an appropriate jumper. 

2.5 Power supply and PCB layout 

As mentioned previously, the intention was to create a high performance anemometer 
design. Part of this overall objective was the desire to minimize electronic noise. For this 
reason the PCB can be configured with dual power supplies. Testing has revealed that 
the broad-band noise level at the TOB is halved when dual supplies are used. Dual power 
supply operation also may be necessary in situations where the bridge requires relatively 
large currents. This could be the case with the optional symmetric bridge configuration 
and/or if probes are used that require a larger current e.g. hot films. Jumpers J1 and J 2 
must be installed if only a single ±15V power supply is to be used. It is recommended 
that the lowest noise modular ±15V power supplies be used. Note that the front panel 
toggle switch labeled POWER does not switch the 110V 60Hz a.c. power. The switch is 
responsible for connection of ±15V power to the PCB. From the viewpoint of minimizing 
noise the advantage of the dual supply configuration is that one supply is dedicated to 
the anemometer section alone (i.e. the bridge, feedback amplifier and power transistors) 
while other ±15V supply is used for the output amplifier where the susceptibility to con- 
tamination by noise is much less. This second ±15V power supply also feeds a single 
+5Y voltage regulator for powering the logic circuits and the front panel indicator LEDs. 
Bypass capacitors are located between the supply and ground rails as close as possible to 
each active device mounted on the PCB. All operational amplifiers use 2 . 2 / / F capacitors 
while 0.1/iF capacitors are used for the digital I.C.s and analog switches. 
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‘Bridge resistors are 3/4 watt metal film 50ppm 
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FIGURE 2.2 Schematic of PCB and front panel. 
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FIGURE 2.2 Schematic of PCB and front panel (continued). 
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3. FABRICATION 


3.1 Chassis and Panels 

Complete assembled bare chassis of the required size are not available commercially but 
the chassis is easily manufactured. The front panel is integral with the PCB since the 
BCD rotary switches are soldered onto the board and the spindles and knobs must align 
with holes in the front panel. 

The top and bottom of the chassis consist of two inch wide aluminum extrusions available 
under part number SR2-0625s-41 from Vector Electronic Company (see Table 1). and must 
be cut to 16 inches length. The edges of the PCB are supported by slots in the extrusions 
as shown in figure 3.2. The extrusions also have holes for fixing the front and rear panels. 

The left and right side panels are cut from 16 gauge aluminum sheet. The right side is 
plain but the left side contains holes for supporting the PCB via standoffs and mounting 
the modular power supplies. The dimensions and location of the holes are shown in figure 
3.6. The side panels engage slots in the top and bottom extrusions and they are free to 
slide, being constrained by only the front and rear panels. Blanks for the front and rear 
panels must be cut to external size from 3/32 inch thick aluminum sheet. After drilling 
the four mounting holes as shown in figure 3.3 the bare chassis can be assembled. 

Figure 3.1 is a drawing of the assembled front panel. The relationship between the 
chassis, PCB and front panel components is shown in figure 3.2. Figure 3.3 shows the 
position of all the holes in the front panel. The location and size of the front panel holes 
is shown in figures 3.4. The specifications for the rear panel holes are shown in figure 3.5. 


3.2 Printed Circuit Board (PCB) 

A complete schematic of the PCB and associated front panel controls is shown in figure 
2.2. The PCB is relatively simple and cheap to manufacture because it has only two 
layers, although it is double-sided. Photocopies of the four photoplot.s used to fabricate 
the PCB are shown in figures 3.7. The specifications for photoplots and construction of 
the PCB are contained in a set of Gerber files which are industry standard files used by 
PCB manufacturers. The following list of files are required to construct the PCB. The files 
are available on a 3.5 inch DOS formatted disk, by sending a written request to: 

Branch Chief 

Fluid Mechanics Research Laboratory Branch 
M/S 260-1 

NASA Ames Research Center 
Moffett Field, CA 94035 


ap_photo.txt 

gpb.cs.ger 

gpb_ss.ger 

pt.cs.ger 

pt_ss.ger 


Text File 
Gerber File 
Gerber File 
Gerber File 
Gerber File 


Files for Photoplot Fabrication 

Aperture for photo plots only 
Component Side Ground Plane 
Solder Side Ground Plane 
Component Side Pads and Traces 
Solder Side Pads and Traces 
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FIGURE 3.1 Layout of assembled front panel with controls labeled. 
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FIGURE 3.2 Relationship between chassis, PCB and front panel components. 
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Drill 0.166 (#19) 
C’sink 0.090 deep 
(4 panel mounting holes) 



FIGURE 3.3 Drawing of front panel holes. External dimensions of panel and position and 
size of holes for mounting to chassis and rack specified. 
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FIGURE 3.4 Position and size of front panel holes for mounting components. 
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FIGURE 3.4 (Continued) Position and size of front panel holes for mounting components. 
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FIGURE 3.5 External dimensions of rear panel and specification for holes. 
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FIGURE 3.7(a) Photocopy (50% reduction) of original photoplot for solder side showing 
traces, pads, ground plane and ground plane block. 

sb.cs.ger Gerber File Component Side Ground Block 

sldr_msk.ger Gerber File Solder Side Ground Block 

ss.cs.ger Gerber File Component Side Silk Screen 

Files for PCB Fabrication 

ap_hole.txt Text File Aperture for PCB Holes Only 

holes. ger Gerber File Hole Only 

Some manufacturers make photoplots for construction of the PCB s from the Gerber 
files while others have this work performed by outside organizations. In either case it is 
cost effective to obtain the photoplots after manufacture so that they can be stored in 
a safe place. Having the photoplots provides a degree of independence from particular 
manufacturers and the cost of having extra boards made in the future is reduced. 
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FIGURE 3.7(b) Photocopy (50% reduction) of original photoplot for component side show- 
ing traces, pads, ground plane and ground plane block. 


The PCB has a silk screen with labels of the components, all pads with holes drilled and 
deburred. a ground plane and solder mask coatings. The cost of fabricating the boards 
was around $60 each for quantity 10. The manufacturer delivered finished boards ready 
for soldering of components. The photoplots for the PCB’s in the sets constructed for 
the FML at NASA Ames Research Center were manufactured independently by Image 
Technology Incorporated. Two sets of PCB’s were manufactured at two different times. A 
list of manufacturers is supplied at the end of this section. 

3.3 Assembly 

Fabtronics Systems was used on two occasions to assemble 20 anemometers and at a later 
date S anemometers, chassis and all components and they delivered complete but untested 
units. This service included plating of the bare chassis and complete preparation of the 
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FIGURE 3.7(c) Photocopy (50% reduction) of original photoplot showing solder mask for 
both sides. 

front panel from the blank provided, the front panel holes and silk screening of the multi- 
colored labels, components were soldered to the PCB and the front panel was attached 
via the BCD rotary switches, and they were installed on the front panel along with other 
hardware such as the BNC connectors, toggle switches, LEDs and potentiometers, and 
power supplies were installed in the chassis and the internal wiring was completed, offered 
by Fabtronics is not unique and should be available at many other locations, chassis for 
mounting in an instrument rack, as shown on page (i). 

3.4 Preliminary Tests and Tuning 

When the sets arrived after fabrication the preliminary test consisted of connecting 
power without a probe attached. If the fuse blows then an ohmmeter should be used to 
check for short circuits between the supply rails. A typical cause of short circuits is poor 
manufacturing techniques in which burrs or other debris cause a short between the power 
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FIGURE 3.7(d) Photocopy (50% reduction) of original photoplot for silk screen labels of 
components. 


supply rails. However this type of error was not found with the 30 PCBs manufactured. 
The most common cause of short circuit was found to be the reversed installation of 
a bipolar capacitor or integrated circuit. When all short circuits are removed the next 
step is to set the magnitude of the square wave to be around ±100 mV by adjusting the 
potentiometer labeled RV 2 located on the PCB. 

It is important that the Gain. Filter (frequency response) and Offset voltage of the 
feedback amplifier are within the range required for stable operation before attempting 
to operate with a hot-wire connected. If the Offset voltage is not within range then the 
anemometer will not operate. If the Filter and Gain settings are too high then the system 
is likely to be unstable. The use of a 1/32 amp fast blow fuse was found to be useful as a 
substitute for a hot-wire for testing purposes. This fuse has a cold resistance of between 
300 and 400 and was operated with a resistance ratio around 2:1. 

Gain setting 3 and filter setting 2 are recommended for this test. Instabilities should 
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be removed by adjusting the Gain, Filter and Offset of the feedback amplifier. If stable 
operation cannot be obtained then each stage of the feedback amplifier should be checked. 
Jumpers J 6 and J7 connect the bridge to the instrumentation amplifier. When these 
jumpers are removed the feedback is disconnected and the amplifier may be tested by 
applying a known input and checking the outputs of each amplifier in the cascade. A 
sinusoidal voltage of 5m V peak-to-peak is recommended because of the high gain of the 
feedback amplifier. 

Similar tests can be performed if problems are experienced with the Output Amplifier. 
Testing and tuning of the Output Amplifier is simpler since the CAL input can be used to 
drive the amplifier directly when the anemometer is in standby mode. 

3.5 Tuning for Maximum Performance 

In order to achieve stable high frequency operation, the feedback amplifier should be 
tuned so that its frequency response is as high as possible. It is also desirable to have 
independent control of the gain and frequency response of the feedback amplifier. However 
the gain and frequency response settings will become related at high frequencies i.e. the 
system will roll-off at higher gain settings irrespective of the filter setting. The feedback 
amplifier should be “opened” for tuning purposes by removing jumpers J6 and J7 which 
connect the bridge to the instrumentation amplifier. A signal generator capable of supply- 
ing a small “clean” signal around, 5mV rms is needed, along with a high frequency rms 
voltmeter. 

Each stage of the amplifier cascade should be tuned in succession by adjusting the capaci- 
tors such that the frequency response is maximally fiat. An HP3325A Synthesizer/Function 
Generator generated the signals used to produce figures 3.8, 3.9 and 3.10. An HP3458A 
Multimeter was used to measure the rms input from the HP3325A and the output of the 
amplifier. Many general purpose multimeters do not a sufficiently high frequency range for 
the rms voltage measurements. The HP3325A can measure rms voltages accurately over 
a frequency range from 1HZ to 10 MHz. 

An example of tuning the feedback amplifier is shown in figures 3.8. Testing was con- 
ducted at 16 discrete frequencies with approximately uniform logarithmic spacing from 
100Hz to 100kHz. The effects of the filters were minimized by using filter setting S for 
these tests. The results for the untuned feedback amplifier are shown in figure 3.8(a) and 
the roll-off can be detected at frequencies as low as 20kHz although -3db point for the 
worst case occurs around 80kHz. This type of untuned behavior of the feedback amplifier 
is acceptable for applications in which the required system frequency response is less than 
30kHz. However tuning is recommended for applications requiring a system frequency 
response in excess of 100kHz. The results of the tuning operation on this particular set 
are shown in figure 3.8(b). The frequency response is essentially flat up to 100kHz. 

The results of tuning another feedback amplifier with filter setting 8 is shown in figure 
3.9. The measurements were performed using 40 discrete frequencies over a wider range 
up to 1MHz. These results indicate that the variable gain stages (the AD744 operational 
amplifiers) are responsible for the lower roll-off frequencies with increasing gain. Assuming 
that each variable gain stage has identical behavior and that the roll-off can be treated 
as a simple pole, then the -3dB point of each variable gain stage can be estimated from 
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FIGURE 3.8 Typical open-loop frequency response characteristics of the feedback amplifier 
(a) Before tuning (b) After tuning. 
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Gain =1 Gain = 4 Gain - 7 

Gain = 2 Gain = 5 Gain = 8 


Gain =1 Gain = 4 Gain - 7 

Gain = 2 Gain = 5 Gain = 8 



Frequency (Hz) 

FIGURE 3.9 Results of tuning operation on another unit showing a wider frequency range. 

the -6dB point of the plots. The amplitude ratio is plotted using a linear scale and the 
-6db points correspond to e 0 /e* ~ 0.5. With each stage set at unity gain, the -6dB point 
(i.e. -3dB for each stage) is in excess of 1MHz. For gain setting 7, each stage has a 
gain A" « 8 and the individual -3dB points occur at a frequency in excess of 200kHz. 
It should be noted that the overall gain of feedback cascade is 6,400 for these settings 
leading to a GBW « 1 x 10 9 Hz. Further, these tests subject the variable gain stages 
to large signals. Analog Devices specifications for the AD744 amplifier is a small signal 
unity gain frequency response of 13 MHz and a full power response of 1.2 MHz. Direct 
testing has not been performed on an isolated stage so the effects of stray capacitance and 
other influences related to the PCB layout have not been determined. However it appears 
that these devices are peforming according to the specifications. It should be noted that 
the fluctuating component of the signal is much smaller when the feedback is enabled by 
connection to the bridge. Therefore the frequency response of the feedback amplifier will 
be higher under operating conditions. 

The effect of the filter controls is shown for gain setting 1 in figure 3.10(a) and gain setting 
7 in figure 3.10(b). The frequency response characteristics of the cascade are essentially 
identical for filter settings of 4 and lower. This is to be expected since filter setting 4 
corresponds to a roll-off frequency of 112kHz which is less than the frequency response of 
the variable gain stages. Minor differences are observed for filter setting 5 corresponding 
to a -3dB point of 225kHz which is comparable to the frequency response of the variable 
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FIGURE 3.10 Effect of filters, (a) Gain setting 1. (b) Gain setting 7. 
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gain stages at gain setting 7 as shown in figure 3.9. The effect of higher filter settings for 
gain setting 7 is much reduced compared to the effect for gain setting 1 since the roll-off 
of the variable gain stages dominate the response. 

3.6 Customization 

The analysis in section 5 suggests that the only way to nullify both the inductive and 
capacitive bridge components is to use a symmetrical bridge. The system may be config- 
ured with an optional symmetrical bridge as described in section 4. The modifications lead 
to R a = R c = 200ft and also require the removal of the Hot-Wire Resistance Thumbwheel 
switch and balance inductor. The purpose of this configuration is to allow the bridge resis- 
tor R f> to be placed at the end of a length of co-axial cable that has been carefully chosen 
to match the probe cable. This particular option has not been tested. 

The gain and filter setting sequence shown in the schematic is approximately logarithmic 
i.e. each sequence involves a doubling operation. This gain and filter sequence ensures 
compatibility with a wide range of probes and operating conditions. However for a given 
probe and operating condition only a small subset of the range of settings can be used 
and it is not possible to obtain optimum performance by fine tuning. If the anemometer 
will only be operated with similar probes over the same range of operating conditions 
then the gain and filter sequence can be modified to suit this range. These modifications 
are accomplished by replacing the input and feedback components of the variable gain 
and filter stages. Using a linear sequence for the gain and filter settings would allow the 
operating point of the anemometer to be more highly optimized. 

The tunable Balance Inductor was selected for use with 5fim diameter wires and a co- 
axial probe cable 5 meters long. The adjustment range is satisfactory for systems with a 
frequency response of 20kHz to 30kHz, but smaller values of the inductor may be required 
for optimum frequency response when smaller wires are used since they have a higher 
frequency response. (Refer to section 5 for suggestions and examples for improving the 
stability when using miniature wires). For optimum performance it is recommended that 
the shortest possible co-axial cable be used for the probe. If a shorter cable is used 
then a smaller Balance Inductor will be required. The scaling for selection of a suitable 
replacement inductor can be based simply on the ratio of the cable lengths. 

The broadband noise level of the TOB output is typically around 2mV to 5m\ rms 
depending on the filter setting of the feedback amplifier. The frequency of the a.c. power 
and its harmonics show up clearly in spectral measurements. Since the anemometer wall 
function with ±12V d.c. supply voltages, tests were performed using automotive batteries 
as a substitute for the T15V modular power supplies. The batteries led to complete 
disappearance of the spikes at the a.c. frequency and its harmonics and reduced the 
broadband noise level by about 10% to 20%. The best signal- to-noise ratio is obtained by 
tuning the frequency response for the required application, and no higher. Unfortunately, 
an increased high frequency response also means an increased susceptibility to noise. 


3.7 Manufacturers, Parts List and Suppliers 

The following list of manufacturers were used for the photoplots, manufacture of the 
PCB’s and fabrication of complete but untested assemblies: 
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Manufacturer of the Photoplots 
Image Technology Inc. 

S21 San Antonio Road 
Palo Alto, CA 94303 
Tel: (415) 494-3113 
Fax: (415) 494-2947 

Manufacturer of second set of 10 PCBs 
Hunter Technology 
1509 Berger Drive 
Santa Claxa, CA 95051 
Tel: (408) 736-5400 
Fax: (408) 736-1908 


Manufacturer of first set of 20 PCBs 
Tel Tec Corporation 
3305 Keifer Road 
San Jose, CA 95112 
Tel: (408) 297-9897 
Fax: (408) 297-1540 

Assembly of 28 Anemometers 
Fabtronics Systems 
4749 Bennett Drive G 
Livermore, CA 94550 
Tel: (510) 447-1316 
Fax: (510) 447-1316 


The following list of suppliers were used for obtaining the components listed in Table 1. 
In some cases the components could only be purchased from the manufacturer. These man- 
ufacturers are Analog Devices, Burr Brown, Elma Electronics, Modular Power Converters 
and PD&E Electronics. Components manufactured by Vector Electronic Company were 
purchased directly from the manufacturer but it may be possible to purchase these items 
through distributors. The other companies listed axe general distributors of electronic 
components. The total cost of these components is under $700. Neither the author nor 
NASA Ames Research Center makes a specific recommendation regarding these suppliers. 
The components are not unique and substitutes are freely available. 


List of Suppliers 

Arrow Electronics 
521 Wedell Drive 
Sunnyvale, CA 94086 
Tel: (408) 745-6600 


Analog Devices 
One Technology Way 
P.O. Box 9106 
Norwood, MA 02062-9106 
Tel: (617) 329-4700 

Burr Brown 
Inti. Airport Industrial Park 
Tucson, AZ 85734 
Tel: (602) 746-1111 

Elma Electronics 
41440 Christy Street 
Fremont, CA 94538 
Tel: (510) 656-3400 

Modular Power Converters Inc. 
RFD Box 41 
Fremont, NH 03044 
Tel: (603) 642-5913 


Digi-Key Corporation 
701 Brooks Avenue South 
P.O. Box 677 

Thief River Falls, MN56701-0677 
Tel: (800) 344-4539 
Hamilton Avnet Electronics 
1175 Bordeaux Drive 
Sunnyvale, CA 94089 
Tel: (408) 920-1211 

Newark Electronics 
1975 Hamilton Avenue 
San Jose, CA 95125 
Tel: (408) 559-6900 
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PD&E Electronics 
Newington Industrial Park 
Newington, NH 03801 
Tel: (603) 436-1211 


Vector Electronic Company 
12460 Gladstone Avenue 
P.O. Box 4336 
Sylmar, CA 91342-4336 
Tel: (818) 365-9661 
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TABLE 3.1 List of Semiconductors and modular electronic components and suppliers. 
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Switches 



TABLE 3.1 (continued) List of switches, connectors and miscellaneous hardware and sup- 
pliers. 
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Total $225.01 





Capacitors 



TABLE 3.1 (continued) List of capacitors and suppliers. 
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Total SI 18.01 








4* OPERATOR GUIDE 


4.1 Overview of Controls 

Figure 3.1 is a schematic of the front panel where labels are used to define the controls 
and LED indicators. The five rotary switches on the left-hand side of the panel have 
BNC coded digital outputs corresponding to ten positions. The digital outputs are used 
to control analog switches on the PCB which connect different resistors across the devices 
which form the adjustable gain and filter blocks. It is important to realize that only the 
first three BCD outputs are used for the filter setting. This means that only the first eight 
positions correspond with the sequential operations described below. (Switch positions 9 
and 10 have the same effect as positions 1 and 2 respectively). 

The uppermost on/off toggle switch and LED indicator for the power refer to the con- 
nection of ±15V power to the PCB. The 110V 60Hz AC power is not switchable and it is 
connected to the power supplies via the rear panel connector and internal fuse. The power 
supplies are responsible for a large proportion of the internal heating so it is preferable 
to leave them on continuously, even when the anemometer is in standby mode, since this 
will maintain the enclosure close to the thermal equilibrium point to be experienced under 
operating conditions. 

The horizontal lines which are inscribed on the panel serve to divide it into five sections. 
The central section contains the controls for the gain, frequency response and offset voltage 
of the feedback amplifier as well as the tunable cable compensation inductor which is 
located in the balance arm of the Wheatstone bridge. The labels on the front panel and 
the buttons for the rotary switches are color-coded red in this region. These controls 
are used to adjust the stability and frequency response of the anemometer as discussed 
in section 2.3. Also located in this section is the Standby/ Operate switch and indicator 
LED. When the switch is in the Standby position the input to the power transistor pair is 
connected to ground. These transistors are responsible for providing current to the bridge 
for heating the wire so that grounding the inputs disables the anemometer. The switch 
must be in the standby position when the probe is connected or disconnected from the 
anemometer otherwise transient currents may cause the probe to burn out. Additionally, 
when the Standby /Operate switch is in the standby position, any signal that is applied to 
the Calibration Input BNC Connector is connected to the input of the output amplifier. 
The Calibration Input BNC connector (labeled CAL) is located in the output amplifier 
section. When the Standby/Operate switch is in the Operate position the output of 
the feedback amplifier is connected to the power transistors and the LED indicator is 
turned on. This connection completes the feedback circuit thereby making the anemometer 
operational. Simultaneously, the Calibration Input connection to the output amplifier is 
replaced with the output voltage of the anemometer. Note that the anemometer output 
voltage is defined to be the voltage at the Top of the Bridge (TOB). The anemometer 
output voltage is available at the BNC connector labeled TOB which is located within the 
group of four BNC connectors on the lower portion of the panel. The connection of the 
TOB voltage to the output amplifier is performed internally on the PCB using an analog 
switch. 
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Immediately below the power switch and LED indicator are the controls for the output 
amplifier section. The labels on the front panel and the buttons for the rotary switches are 
color-coded blue in this region. Typically, the output voltage of a constant-temperature 
hot-wire anemometer consists of a relatively large d.c. component upon which are super- 
imposed the small fluctuations which are of interest. As mentioned in section 2.4, it is 
common practice to utilize a separate “buck and gain 1 ’ amplifier to subtract a d.c. volt- 
age from the anemometer output and to amplify the resultant signal so that the range of 
voltages fills the voltage window of a digitizer. For example, the numerical resolution of 
analog-to-digital converters are limited and the output amplifier can be used to maximise 
the numerical resolution of the fluctuations. 

For convenience, an output amplifier has been incorporated into the instrument for this 
purpose. The amplifier is essentially an analog computer circuit which is hardwired 
to perform the linear transformation given by equation (2.1). The output amplifier is 
configured as an inverting amplifier because the sensitivity of the anemometer (TOB) 
output voltage is negative. With this arrangement the sensitivity of the output voltage is 
positive with respect to velocity. The gain G is adjusted using the rotary switch labeled 
GAIN and the gain follows the logarithmic sequence 1, 2, 4, 8 ... 256 which correspond 
to switch positions 1 to 9. Gain setting 10 also has a gain of 256. The offset voltage E 0 n 
is adjusted using the locking counting dial labeled OFFSET on the panel. As mentioned 
above, when the system is switched to Standby mode the input voltage Ei consists of the 
signal applied to the Calibration Input BNC connector. This allows the constants G and 
Eoff in equation (2.1) to be determined by direct calibration. One reason for wanting these 
constants with high precision is to enable the recalculation in software of the TOB voltages 
from the digitized values of the output amplifier signals. This is necessary for hot-wire 
calibration schemes that are based on King’s Law for example. The output amplifier is 
configured with a low-pass filter consisting of two simple poles in series which are controlled 
by the rotary switch labeled FILTER. The -3dB points follow the logarithmic sequence 
shown on the schematic in figure 2.2. It should be emphasized that these filter values are 
nominal and that tuning capacitors are provided for adjustment of the roll-off frequency 
and to ensure that the Bode diagram is flat. The Bode diagram cannot be guaranteed to 
be flat out to these -3dB points unless the tuning capacitors have been tuned manually. 
The one exception corresponds to switch position 8 since the tuning capacitors have the 
smallest effect in this position. However tuning capacitors are provided for setting the 
maximum frequency response and these should be adjusted if the signal content exceeds 
around 25kHz. It is highly recommended that the Filter switch be used in position 8 unless 
the circuit has been manually tuned for the other filter settings. 

No provision has been included for measuring the cold resistance of wires because this 
function should only need to be performed occasionally. Inclusion of the extra electronics 
for such an infrequently used function was considered to be of the lowest priority early in 
the design stage. Finally, it became apparent that there would be insufficient space for 
the extra controls despite the rather large size of the front panel. Therefore some other 
device must be used for measuring the cold resistance. This device should be selected with 
care. Some Ohm meters will pass a small but significant current through the wire and of 
course this will change the resistance of the wire. A Kelvin bridge or other null current 
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device is best for this purpose. The operating resistance of the hot-wire is set using the 
Hot-Wire Resistance Thumbwheel Switch labeled RESISTANCE in figure 3.1. The switch 
is configured such that the setting corresponds to the balance resistor setting in ohms i.e. 
the setting shown in figure 3.1 corresponds to R(, = 123.0f2. In the standard configuration 
the cross-bridge ratio is 10:1 so that this thumbwheel switch setting corresponds to a hot 
operating resistance of R u ,=l2.3fl for the hot-wire. 

A square- wave test is required for tuning the anemometer and this capability is provided 
via the Square Wave toggle switch. The LED indicates when the square- wave circuit is 
turned on. In order to minimize noise, all the digital circuitry used for the square- wave is 
powered down when the switch is the off position. The amplitude of the square- wave signal 
must be adjusted on the PCB using the trimming potentiometer RV2. The amplitude 
should be small enough to avoid nonlinear behavior i.e. the response to postive and 
negative fluctuations should be approximately symmetrical. The square-wave frequency is 
not a critical parameter and the rotary control switch is provided merely for convenience. 

4.2 Operation 

The hot-wire should be connected to the BNC connector labeled PROBE which is located 
among the group of four connectors on the lower portion of the front panel as shown 
in figure 3.1. As mentioned previously, the anemometer output voltage TOB and the 
output amplifier voltage are also available in this region. If the anemometer appears to 
be inoperable then the use of 1/32 amp fuse is recommended for testing purposes. Use 
of the fuse may avoid burning out a good hot-wire. Section 3.4 “Preliminary Tests and 
Tuning” describes some of the procedures used to set the Offset Voltage and Gain and Filter 
settings. Experience has shown that the most common cause of problems is associated with 
an inappropriate Offset Voltage setting for the feedback amplifier. 

After the anemometer is operational the next step involves setting the stability and 
frequency response. The toggle switch in the region labeled SQUARE WAVE on the front 
panel introduces a square wave perturbation to the Offset Voltage of the feedback amplifier. 
Tuning of the anemometer should be performed by observing the square wave response 
at the TOB output on an oscilloscope as the controls are adjusted. Using the output 
amplifier for tuning purposes is not recommended because of the extra time constants 
associated with this amplifier. Optimal setting of the square- wave response is discussed in 
section 6. It is emphasized here again that the electronic square-wave test is an indirect 
method for setting the stability and frequency response. Proper interpretation of the 
square- wave response often requires subtle insights into the system behavior. It is easy for 
an inexperienced operator to be deluded and infer a frequency response that is much higher 
than the actual system response. Examples are provided in Section 6 which demonstrate 
how this can occur. Additional cases are considered that should provide a useful guide for 
relating the influence of the anemometer controls on the square-wave response. The advice 
contained in this section should also prove useful to experienced operators. 

The BNC connector labeled EXT RB is required for the optional configuration in which 
the cross-bridge ratio is altered from 10:1 to a 1:1 symmetrical bridge by setting jumpers 
on the PCB. With this optional symmetrical bridge configuration the Hot-Wire Resistance 
Thumbwheel switch is removed from the bridge and the balance resistor Rb must be located 
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at the end of a length of co-axial cable that has been carefully chosen to match the probe 
cable. This configuration offers superior compensation to the simple lumped inductor 
(by using actual cable instead). Once again, it is emphasized that the bridge must be 
symmetrical for this option to work properly. 

The Cable Compensation Inductor was selected for use with 5/im diameter wires and 
a 5 meter long probe cable. The adjustment range is satisfactory for systems with a 
frequency response of 20kHz to 30kHz, but smaller values of the inductor may be required 
for optimum frequency response when smaller wires are used since they have a higher 
frequency response (see section 6.5). One simple solution in this situation is to use a 
longer cable rather than substitute a smaller inductor. However the overall performance 
of the system will be improved if a smaller inductor is used. In this case the benefits of 
reducing the cable length to a practical minimum should be considered in conjunction with 
the substitution of the inductor since performance is improved when the total inductance 
is reduced. The scaling for selection of a suitable replacement inductor can be based on 
the ratio of the cable lengths. 
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5. ANALYSIS 


5.1 Introduction 

In the constant-temperature mode of operation the hot-wire filament is located in a 
Wheatstone bridge and a feedback amplifier is used to maintain the wire at very nearly 
constant resistance. It will be shown by example in section 6 that the feedback amplifier 
must have the potential for both a high gain and a high frequency response in order for 
the system as a whole to achieve stable high frequency operation. Practical constant 
temperature hot-wire anemometer designs should use a cascade of amplifiers rather than 
a single feedback amplifier because of the gain- bandwidth product {GBW ) limitations 
of physical devices. Therefore a cascade of n amplifiers is used in the model shown in 
figure 5.1. The device labeled I is a current booster stage which usually consists of a pair 
of power transistors. It will be assumed that the current booster has unity gain and zero 
offset exactly. The assumption of a perfect voltage follower means that the current booster 
need not be considered in the analysis. The behavior of each amplifier will be assumed 
to be linear with zero offset voltage. The net offset voltage of the amplifier cascade is 
controlled by introducing a d.c. voltage E qi and a perturbation voltage es into the 777 th 
amplifier. The d.c. voltage is used to adjust the frequency response and stability of the 
system while the perturbation voltage is used for estimating the system frequency response. 
The perturbation voltage usually consists of a square- wave although sinusoidal inputs are 
sometimes used. In general, upper case letters refer to d.c. quantities while the lower case 
letters refer to small perturbations. Z a (s ), Zb{s), Z c (s) and Z w (s) represent the electrical 
impedances of each arm of the bridge and £ is the Laplace variable. 

5.2 The static operating point 

Only the resistive components of the Wheatstone Bridge need to be considered for the 
static analysis below i.e. Z a (s) = i? a , Zb($) = Rb , Z c (s) = R c and Z w (s) = R w . The 
voltage at the top of the bridge i.e. the output of the n th amplifier is given by 

Eo n = Kb [K a (IlRa — I 2 Rc) + Eqi] (5-1) 

where K a = Ky...K m -i is the total gain of the cascade up to (but not including) the m th 
amplifier where the offset voltage is injected and Kb = K m ...K n is the total gain of the 
cascade from the m lh to the 77 th amplifier. The bridge voltage can also be expressed in 
terms of the currents I\ and J 2 , 


E 0n — I\{R a + Rw) 

(5.2) 

E 0n = l2{Rb + Rc) 

(5.3) 


Substituting E 0n from equation (5.2) and J 2 from equation (5.3) into equation (5.1) leads 
to the following expression for the wire current 

j KbE qi {Rb + R c ) ^ ~ ^ j 

1 ” (R a + R w )(R b + R C ) + KR 
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FIGURE 5.1 Model of Constant-Temperature Hot-wire Anemometer. 

o 

where R = R U ,R C — R a Rb represents the bridge imbalance and K = I\ a Kb is the overall 
gain of the cascade. Equation (5.4) represents the behavior of the circuit. Following 
Smits, Perry and Hoffman (1978), an equation which represents the balance between the 
heat transfer and the Joule heating of the filament is given by 


( \ 1 \ nlK g Nu 

\ Rw J c 


(5.5) 


where R g is the wire resistance at fluid temperature, l is the length of the wire, K g is the 
thermal conductivity of the fluid, Nu is the Nusselt number and c is the temperature coef- 
ficient of resistivity of the wire material. The operating point is defined by the intersection 
of equations (5.4) with (5.5) which must be solved by iteration. 

A number of important asymptotes and limits exist with respect to the circuit equation. 
For example, the r.h.s. of equation (5.4) must be indeterminate if I\ is to remain finite as 
E g i -* 0 i.e. 

( Rq 4“ Rwa ){R b + R c ) + KR a =0 (5.6) 

O o 

where R a = R wa R c — R a Rb is the value of R corresponding to the asymptotic value of the 
wire resistance 

D _ ER a R b — R a {Rb + Rc ) /r 

U ’ a “ KR C + (Rb + R7) (5 '° 

It is important to realize that perfect bridge balance can be achieved only in the limit 
as the offset voltage E qi — > 0 and as the overall gain of the amplifier cascade K —> oc. 

c 

However the quantity KR always remains finite, even as the limits mentioned above are 
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approached. The quantity I\R has a strong effect on the system stability and frequency 
response and it is controlled by adjusting the offset voltage. 

Wood (1975) performed an analysis of the type of system in which the gain of the 
feedback amplifier is frequency dependent. For systems with a frequency dependent gain 
the term I\ R is mult iplied by the ratio of the a.c. gain to the d.c. gain which can be as small 
as 0.003. This minimizes the effect of offset voltage and the system stability and frequency 
response are adjusted by varying the a.c. gain in conjunction with the balance inductance 
and roll-off frequency of the amplifier. This is typical of some commercially available 
systems and is thought to offer several advantages over the type of system considered 
in this report. Firstly, tuning the system frequency response by adjusting the a.c. gain 
does not alter the hot-wire calibration as does offset voltage adjustment in the system 
described here. However if the frequency response and stability of the system are properly 
tuned then there should be no need for further adjustments after calibration. Secondly, 
as mentioned above, a constant-temperature hot-wire anemometer cannot maintain the 
bridge in perfect balance. The wire resistance varies slightly with mean velocity so that 
the calibration curve deviates from the ideal constant temperature behavior. Systems with 
high d.c. gain amplifiers are designed to minimize the bridge imbalance so that the hot- 
wire system can be thought to more closely follow some idealized heat transfer law, such as 
King’s law. Smits and Perry (1980) studied the influence of the overall gain by assuming 
King’s law and solving equations (5.4) and (5.5) for the operating point as the velocity is 
varied. They found that the deviations from ideal behavior were of order 1% for a gain 
of 500 and that increasing the gain beyond 1000 brings only a diminishing return. These 
sm all deviations from ideal behavior would be accounted for in practice by the curve fit 
used for the calibration. In any case, a hot-wire filament cannot be regarded as a straight 
infinitely long cylinder of uniform circumferential and spanwise temperature (see Perry 
1982) so it seems unlikely that there exists a single universal heat-transfer law for hot-wire 
calibrations. W hen these practical considerations are taken into account the advantage of 
frequency dependent gain systems should be considered idealistic and of minor significance. 

Of more concern is that a negative step in the Bode diagram occurs when the feedback 
amplifier has a frequency dependent gain. One shortcoming of the analysis of Wood (1975) 
is that it applies only for frequencies below which the a.c. gain is flat. Although he could 
not calculate the response from his analysis in the frequency range where the amplifier 
gain is falling from its d.c. level to its a.c. plateau, Wood realized that a significant step 
in the frequency response could occur. He argued, by worst case example, that the effect 
is only important (7%) at low resistance ratios and for low a.c. gains. However Smits 
and Perry (1980) found that the step size can be significant for resistance ratios greater 
than 1.1. Even for resistance ratios of order 1.5 they showed that the size of the step can 
be as large as 3% of the static response. For these reasons, the feedback amplifier of the 
anemometer described in this report was designed to have a flat frequency response. 

5.3 Transfer functions for dynamic response 

The circuit shown in figure 5.1 will be analyzed using Laplace transform methods since 
transient phenomena can be examined (such as the square-wave response) in addition to 
the usual steady-state sinusoidal response. It will be assumed that the dynamic behavior 
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(5.8) 


of any j th amplifier in the cascade can be adequately described by 



Kj 


Ms) 

Bj(s ) 


which is the transfer function for output voltage fluctuations in terms of the input voltage 
fluctuations. The circumflex denotes the Laplace transform i.e. 




r 00 

= L{e;.(/)}= / ei-(t)e si dt 

Jo 


(5.9) 


Kj is the amplifier gain and the roots of the polynomials Aj(s) and Bj(s) are the zeros 
and poles of the amplifier. An overview of the application of Laplace transform methods to 
systems of differential equations is given in Appendix A along with definitions of a transfer 
function and its poles and zeros. From now on the circumflex will be dropped but the 
Laplace transform is still implied. 

In a physical circuit the current booster stage usually has a frequency response well in 
excess of the amplifiers. For the purpose of modeling it will be assumed that the frequency 
response is infinite in addition to the assumptions of exact unity gain and zero offset 
voltage. As mentioned previously, the assumption of an “ideal” current booster stage 
means that it need not be considered in the following analysis. 

The small perturbation output voltage of the first amplifier in the cascade is given by 


j- M j) 


iiZ a (s) - i 2 Z c (s ) 


(5.10) 


where i\ and i 2 are the small perturbation bridge currents and Z a (s) and Z c (s) repre- 
sent the electrical impedances of the upper half of the bridge (see figure 5.1). The small 
perturbation bridge currents ii and i 2 are dependent on both offset voltage and velocity 
perturbations. The system transfer functions will be derived for the two cases of offset 
voltage perturbations (with no velocity fluctuations) and velocity fluctuations (with no 
offset voltage perturbations). 


5.3.1 System Transfer Function for Offset Voltage Perturbations 


In the absence of velocity fluctuations, the small perturbation currents i\ and i 2 can be 
expressed in terms of the voltage fluctuation at the top of the bridge i.e. the output e Qn 
of the last (n th ) amplifier 


Z a {s) + Z w (s) 


(5.11a) 


and 


*2 


Zb(s ) + Z c (s ) 


(5.116) 


The electrical 
i.e. 


impedance of each arm of the bridge can be expressed as transfer functions 


Z a (s) 

Zc(s) 


ZaM 
Z>a D {s) 
ZcA$) 
Z Cd (sY 


and 


Zb($) 

Z w {s) 


ZbM_ 
Zb D {s)' 
Z WN (^) 

Z WD (s) 


(5.12) 
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This representation is completely general and an arbitrarily large number of elements can 
be used. The only restriction is that the elements must consist of lumped components. 
Using equations (5.11a) and (5.11b) to substitute for and e? in equation (5.10) and then 
applying equation (5.S) along the amplifier cascade leads to the expression for the output 
of the last (n th ) amplifier <„ n . i.e. 


e 


On 


-4m( , s)'''-4. n (.s) 
" B m (s) ■ ■ ■ B n {s) 


Ki- 


■ ■ K m -t 


.4i(s) • • • 

B](S ) • • • B,r,-i(s) 


Z a (s) + Z w (s) 


c(s) + Z b (s)J 


(5.13) 


Equation (5.13) implies that a minimum of two equivalent amplifier stages are needed to 
properly account for the introduction of offset voltage perturbations. The first of these 
two equivalent amplifiers, a, has a transfer function equal to the amplifier cascade up to 
(but not including) the 777 th amplifier where the offset voltage is injected, i.e. 



e i« 


• • • K m -\ 


4.i(s) • • • ^4 m -i(s) 
B\{s) ■ ■ • £? m _i(s) 


= K 


a 


Ms) 

B a (s ) 


(5.14) 

(5.15) 


The second equivalent amplifier, 6, has a transfer function equivalent to the amplifier 
cascade from the m th to the n th amplifier, i.e. 



T r ^m(^) * * ' 

n TZ 7 7 ZT 7 7 


= K b 


4.»(s) 

B b (s) 


(5.16) 

(5.17) 


The transfer function of the system response to offset voltage fluctuations using the equiv- 
alent amplifiers a and b is obtained from equation (5.13) as, 

77 w = a 'pw (5 ' 18) 

where K e = I\ b and the polynomial expression for the system zeros is given by 

Qe(s) = A b (s)B a (s)Z 1 (s) (5.19) 


and the polynomial expression for the system poles is given by 

P(s) = B a {s)B b (s)Zi(s ) + K a K b A a (s)A b (s)Z 2 (s) 


(5.20) 
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FIGURE 5.2 General representation of Bridge Impedance Z w required used to derive the 
System Transfer Function for Velocity Fluctuations. 

where. 

Zi(s)=[Z aN (s)Z WD (s) + Z aD (s)Z WN (s)][Z bN (s)Z CD (s) + Z bD (s)Z CN ( S )} (5.21) 

and 

^(•s) = ( s )Zt, D (s)Z CN (s)Z WN (s) — Z aN (s)Zb N (s)Z CD (s)Z WD (s) (5.22) 

Note that the poles of the amplifiers which precede the offset voltage injection stage appear 
as zeros in the overall system transfer for offset voltage fluctuations in equation (5.19). 

5.3.2 System Transfer Function for Velocity Fluctuations 

Voltage fluctuations across the hot-wire filament originate from both velocity fluctuations 
and perturbations in wire current. In the previous case, the fluctuations in wire current 
were the result of offset voltage perturbations and the current perturbations in each arm 
of the bridge could be simply expressed in terms of the bridge impedances in equations 
(5.11). The voltage fluctuations at the top of the bridge could then be calculated directly 
using equation (5.13). In the case of velocity fluctuations the wire acts as the source of the 
voltage fluctuations and the current perturbations depend on the nature of the impedances 
forming the wire arm of the bridge Z w (s). The representation of Z w (s) required for the 
derivation of the transfer functions for velocity fluctuations is shown in figure 5.2. The 
operational impedance Z s (s) is in series with the wire and Z p (s) is the impedance in 
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parallel with the wire and Z s (s). This representation is also completely general and an 
arbitrarily large number of elements can be used to form the impedances Z s [s), Z w j re (s) 
and Z p (s). 

Let the transfer functions representing the impedances be given by. 


Z w \ie 


00 = 


Z wire jv ( ^ ) 
Z wire d ( ^ ) 


Z*(s) 


Zs N ( s) 

Z 3D (^) 


and Z p {s) 


Z ps (s) 

z PD (*y 


(5.23) 


The voltage perturbation across the wire is given by 


Cwire — ^wire(^) ?1 a T S u F u [s) U 


(5.24) 


where i\ a is the wire current, S u is the sensitivity to velocity perturbations (i.e. dE Wire / dl ) 
at zero frequency, u is the velocity perturbation and F u (s) represents the normalized 
frequency response characteristics of the wire i.e. F u ( 0) = 1. 5 U is a constant at fixed 
velocity. The voltage across the wire arm of the bridge is given by 


e i 


Zs{$) + Z w i re (s) 


l\ a T Su Xi 


(5.25) 


which is also equal to the voltage across the parallel impedance i.e. 


6\ — Zp(s) i\b . 


(5.26) 


The total perturbation current (ij = i\ a + iib ) can be found from equations (5.25) and 
(5.26) so that the voltage perturbation at the top of the bridge can be determined i.e. 


e 0 „ = 


Z a (s) T Z w (s) 


i\ + S u F u (s)F z (s) u 


(5.27) 


where 

F (s)= Zp{s) 

H ' Z p (s) + Z.(s) + Z wiTe (s) 

and the impedance of the wire arm of the bridge is given by 

„ , x _ %p{ s )[Zs{s) ~t~ ■^wire('S)] 

Z p (s) + Z s (s) + Z wire (s) 


(5.28) 


(5.29) 


Substituting for from equation (5.27) and i 2 from equation (5.11b) into equation (5.10) 
and using the two amplifier model for the cascade, the system transfer function for velocity 
fluctuations is given by 

—(*) = *'. (5 ' 30) 

u *( s ) 


where 


K u = S u K a K b 


(5.31) 
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and 


Qu(s) — -4 a ( s ).4(, ( s )Z :i (s)F u ( s)Z w \ Tt ( s ) 


(5.32) 


where 


■^3 ( ** ) — Z aN ($)Z PN {s)Zft D (s) 


Zbsi S )ZcD ( S ) + -Zfc D (s)Z Cv (,s) 


(5.33) 


The behavior of a hot-wire filament is very complex. For example. Perry ( 1982) discusses 
aeroelastic effects and perturbations in the symmetry of the temperature distribution along 
the wire. Even relatively simple models which allow for a non-uniform temperature distri- 
bution along the wire must use a partial differential equation in space and time. In order 
to keep the analysis tractable a simple lumped model will be used for the hot-wire filament 
i.e. 


and 


Fu{s) = 


1 

T w s + 1 


^wireiy(^) R W T W S + R w -(- Q 

•Zwire£>('S) T W S T 1 


(5.34) 

(5.35) 


T w is the lumped time constant of the wire arising from its thermal inertia, R w is the 
d.c. wire resistance and a = R W (R W — R g )/R g where R g is the wire resistance at fluid 
temperature. The simple lumped model for the filament possesses the same time constant 
for the simple poles of the sensitivity to velocity and current fluctuations and this leads to 
a simplified expression for the system zeros 


Qu{s) = Aa{*)At(s)Z 3 (s) (5.36) 

In the limit Z p (s) — > oo (i.e. when there are no elements in parallel with the wire) 
the expression for Z 3 (s) given by equation (5.33) should not be used since it implies that 
Q„(s) — ► oo. The expressions for Zi(s) and Z 2 (s) given by equations (5.21) and (5.22) can 
be factored by Z PN (s) and when Z p ( $ ) ~ oo the expression for Zj ( ,s ) is given by 


Zz(s) — Z aN (s)Z SD (s) 


Zf, N (s)Z CD (s) -f Zf, D (s)Z Cy (s) 


(5.37) 


Equations (5.21) and (5.22) can still be used to evaluate Zi{s) and Z 2 {s) when Z p (s) = oo. 

It is important to note that the poles of the transfer functions for both offset and velocity 
fluctuations are identical (provided that F u (s)Z w i reD (5) = 1). This result is completely 
general in the sense that the system poles are identical, independent of the number of 
amplifiers in the cascade, the complexity of the overall transfer function for the feedback 
amplifier and the number of lumped components used to model the bridge impedances. 

5.4 Some illustrative examples of transfer functions 

Unlike the poles, the zeros of the transfer function for offset and velocity fluctuations are 
quite different. It has already been noted that the poles of amplifiers preceding the offset 
voltage control stage appear as zeros in the transfer function for offset voltage perturbations 
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but not for the velocity fluctuations. Other differences between the zeros of the two 
system transfer functions depend on the nature of the bridge impedances and these will 
be considered in the examples below. It should be emphasized that the derivation of the 
polynomial expressions for the system poles and zeros in equations (5.19), (5.20) and (5.36) 
can be quite involved, even for relatively simple configurations. It is instructive to consider 
some particular cases. 


5.4-1 The Simplest Possible Configuration 

For the simplest possible configuration only resistive elements appear in the bridge i.e. 


Z a *(s) = R a , 

Z aD ( s ) ~ 1 

(5.38) 

Z bN (s) = R b , 

Z bD (s) = l 

(5.39) 

Z CN (s) = R c , 

Z CD (s) = l 

(5.40) 


and from equations (5.35) and (5.12) the impedance of the wire arm of the bridge is given 

bv 


Z ws (s) — R u ,TxiiS + “b ol 

Z WD (s) = T u ,s + 1 (5.41) 

Substituting equations (5.38) ... (5.41) into (5.21), (5.22) and (5.37) leads to the expres- 
sions 

Zi($) = ( R b + R c )[(Ra + R W )T W S + (Ra + Rw + «)] (5.42) 

^ 2 ( 3 ) = R(T w s + 1)4* R c <y (5.43) 

Z Z (s) = Ra(Rb + Rc) ( 5 - 44 ) 

Assuming an ideal amplifier i.e. one with infinite frequency response, then 

A a (s) = B a (s) = A b {s) = B b (s ) = 1 (5.45) 

Let the feedback amplifier consist of a single stage such that A„ = 1 and K b = A. 
Substitution of equations (5.42) ... (5.45) into (5.19), (5.20) and (5.36) leads to 

Q t (s) = (Ra + R w )(Rh + Rc)T w s + (R a +R W + a)(R h + R c ) (5.46) 

Qu(s) = R b + Rc ( 5 - 47 ) 

P(s) = [(R b +R c )(R a + R w ) + KR]T w s 

+ (R b + R c )(R a + R w + a) + K(k + R c + a) (5.48) 


The time constant of the simple zero in equation (5.46) is approximately the same as the 
time constant T w of the wire filament. It is this zero which causes the system response to 
a step input of offset voltage to appear like the response expected from a delta-function 
impulse of velocity. Note that this system possesses only one simple pole so that the step 
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response can only exhibit an exponential decay i.e. the system response cannot exhibit 
ringing. The static response and the dynamic response are coupled through the quantity 
R which is controlled by adjusting the offset voltage. Consideration of equations (5.6) 
and (5.48) indicates that the time constant of the simple pole approaches zero (i.e. the 
frequency response of the system becomes infinite) as the offset voltage approaches zero. In 
practice the frequency response is limited by other higher order effects. It was once believed 
that the frequency response characteristics of the feedback amplifier were the dominant 
factor in the determination of overall system performance and stability. However Davis 
and Davies (1968) demonstrated the importance of the inductance of the probe cable and 
the offset voltage. Perry and Morrison (1971) extended this work and the results for an 
inductive system are given in the next section. 

5-4-2 Hot-wire System, with Inductance 

Perry and Morrison modeled the hot-wire cable as a simple lumped inductive element 
L u . in series with the wire. Equations (5.35) and (5.12) give the components of impedance 
of the wire arm of the bridge as, 

Z WN (s) = L w T w s 2 -fi {R W T W -fi L w )s -fi ( R w -fi ») 

Z WD (s) = T w s + 1 (5.49) 

They demonstrated that the system stability and frequency response could be improved 
by incorporating a tunable balance inductor Lb in series with the bridge resistor Rb- From 
equation (5.12), the components of the impedance of this arm of the bridge are, 

Zb N (5) = LbS -fi Rb 

Zb D (s) = 1 (5.50) 

The upper bridge elements are purely resistive and substitution of equations (5.38), (5.40), 
(5.49) and (5.50) into equations (5.21), (5.22) and (5.37) leads to the expressions 

Z\{s) -- LbL w T w s 3 

+ { [{Rb + R C )L W + (R a -fi R w )Lb]T w -fi LbL w }s 2 
+ {{Rb + 72 c )[(i?a + Rw)T w -fi L w { -fi ( R a -fi R w -fi a)Lb}s 
+ {Rb + R c ){R a + R w + a) (5.51) 

Z 2 (s) = (R C L W - R a L b )T w s 2 + {RT W + R C L W - R a L b )s + {R + R c a) (5.52) 
Z^{s) = R a [its + (Rb + R c )] (5.53) 

Perry and Morrison did not consider the frequency response characteristics of the feedback 
amplifier so that the poles and zeros are given by equation (5.45). Substitution of these 
expressions into equation (5.20) for P(s) leads to a 3 rt, -order polynomial for the system 
poles 

P(s) = P 3 s 3 + P 7 s 2 + P lS + P 0 (5.54) 


52 



where the constants Pq ...P 3 are given by 
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/Co 
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0 

0 

Cl 

l 0 

0 

1 

0 

Co 

Vo 

0 

0 

\) 

\Cz 


n 0 
° 1 
I 0 0 
Vo 0 



( C'h 0 
Ca-i 
Ca 2 


(5.55) 


Co • • • C 3 and Cao • • • C'a’o are constants that depend on the system parameters i.e. 


Co — (Rb + R c ){Ra + Rw + o ) 

Ci = ( Rb + Rc ) [(f?a + Rtv )T W + £„.] + (R a + Ru- + a )Lb 
C2 = [(-Rt + Rc)Lw + (Ra + Rw )Lb\ T w + LbL w 
C3 = LbL w T w 
C ao — R + Rc<* 

Cai = RT W + R C L W — R a Lb 
C K2 = (RcL w — R a Lb)T w 


(5.56) 


The presentation of equation (5.55) using matrix notation has the advantage of being 
compact while retaining the overall appearance of equation (5.20) and allows the important 
features of the coefficients to be examined in equation (5.56). 

A single zero occurs in the transfer function for velocity fluctuations i.e. 

Qu( s ) = Ra(Lb s + Rb + Rc)- (5.57) 


In most configurations this zero is located well beyond the frequency range of interest so 
that it is not of great significance. The expression for the zeros of the transfer function for 
offset voltage fluctuations is 3 rd -order 

Q e (s) = C iS 3 + C 2 s 2 + C lS + C 0 . (5.58) 

The simple zero given by equation (5.46) for the simple system with only resistive bridge 
elements has a time constant that is approximately the same as the wire filament. It turns 
out there is always a simple zero in equation (5.58) with a time constant of similar value. 
Therefore the response to a step input will also appear like the response expected from a 
delta-function impulse of velocity. However the expression for the poles is now 3 rd -order 
so that a pair of complex conjugate poles are possible. In this case the step response may 
exhibit ringing i.e. the decaying sinusoidal response which users of hot-wire anemometers 
are familiar with. 


5. 4-3 Effect of the Frequency Response of the Feedback Amplifier 

The positions of the system poles on the s-plane are independent of how the poles and 
zero are distributed among the first and second equivalent amplifiers. However extra zeros 
will appear in the transfer function for offset voltage fluctuations if the poles are located 
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FIGURE 5.3 Sketch of Bode diagram for amplitude of a feedback amplifier with a simple 
zero with time constant M\ and two simple poles with time constants M-> and M 3 


before the offset voltage control stage i.e. B a (s). see equation (5.19). The number of zeros 
in the transfer function for velocity fluctuations remains unaltered in this situation, see 
equation (5.36). Further, a zero that is incorporated into the feedback amplifier will only 
appear in both system transfer functions if it is added to the second equivalent amplifier. If 
any poles and/or zeros are deliberately added to the feedback amplifier for control purposes 
then it is highly recommended that they be added only to the second equivalent amplifier 
i.e. after the offset voltage injection stage. This will avoid introducing any further dif- 
ferences between the system transfer functions for velocity and offset voltage fluctuations. 
The inherent differences between the transfer functions already make it difficult enough to 
correctly infer the system frequency response from square- wave tests. 

The addition of a simple zero with time constant M\ and two simple poles with time 
constants A/2 and A/3 to the second equivalent amplifier (equation 5.17) leads to a system 
in which the polynomial for the poles is 5 th -order i.e. 


P(s) = P5S 5 + P 4 s 4 + P 3 s 3 + P 2 s 2 + P lS + Po 


(5.59) 
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where the constants Pq.. .P5 are given by 
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The matrix notation in equation (5.60) illustrates the coupling between the gain and the 
poles and zeros of the feedback amplifier. The expressions for the zeros can be obtained 
by multiplying the polynomials in equations (5.57) and (5.58) by (Mis + 1). As a matter 
of interest this configuration could be used for examination of systems in which the gain 
of the feedback amplifier is frequency dependent by putting M 2 > M\ > M 3 . The Bode 
diagram of a feedback amplifier with this particular frequency dependent gain distribution 
is sketched in figure 5.3. Smits and Perry (1980) considered the behavior of systems with 
a feedback amplifier with a frequency dependent gain but did not account for the high 
frequency roll-off i.e. M 3 = 0. For the reasons mentioned previously this type of system 
has not been used for the anemometer design described in this memorandum. 



(5.60) 


5-4-4 Effect of additional reactive bridge components 

Including the frequency response characteristics of the feedback amplifier introduces a 
relatively small increase in the effort required to expand the transfer function polynmials. 
However the addition of a small number of reactive bridge elements leads to a much higher 
level of complexity. Using a single lumped inductor is a rather simplistic representation of a 
co-axial hot-wire cable with its distributed inductance, capacitance, resistance and leakage 
conductance. A slightly more sophisticated model for the cable could include a lumped 
capacitor C w in parallel with lumped inductor L u . and the wire. A balance capacitor 
Cb could be introduced in parallel with the balance resistor and balance inductor. This 
configuration was studied by Watmuff (1987) and leads to a system in which the polynomial 
for the poles is 7 th -order i.e. 


P(s) = P 7 s 7 + P 6 s 6 + P 5 s 5 + P 4 s 4 + P 3 s 3 + P 2 s 2 + Pis + P 0 (5.61) 
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where the constants Pq...P- are given by 
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The constants Co, C\, C 2 , C 3 , Ca'o, Cai aJid Ca ’2 are given in equation (5.56) and 


(5.62) 


Ci = R b R c (R a + Rw + ot)C b 

+ R a (R b + R C )(R W + ck)C w 
C 2 = R a {Rb ■+■ Rc^RijuTw + L W )C W 
+ R a (R w ■+• Oe)L b C w 
+ R b Rc[{Ra + R W )T W + L w ^C b 
+ R c (R a + Rw + <x)L b C b 
+ R a R b R c (R u , + a)C b C w 
C 3 ' = R a (R b + R c )L w T w C w 
+ Ra{RwT w + L w )L b C w 
+ R c [( R a + R w )Tw + Tu>] L b C b 
+ R b R c L u ,T w C b 


+ R a RbR c (RwTw + L w )C b C w 
+ R a R c {Rw + a)L b C b C w 
C4 = R a Rc [(-RtoTj + R b L w )T w + 

+ ( R c C b + R a C w )L b L u ,T w 
C b = R a R c LbL w T w C b C w 
Chi' = R b (R w + a)(R c Cb — R a C w ) 

Ca'2 = [L b {R w + a ) + Rb(R w T w + L w )\(R c C b — R a C u ) 
CaV = [(RwL b + RbL w )T w + L b L w ](R c C b - R a C w ) 
Cka = L b L w T w (RcC b — R„C W ) 


(5.63) 



The zeros for offset voltage perturbations are given by 


Qe(s) = (Mjs + 1 ) [ C 5 ’s 5 + CU’s 4 + (C 3 + C 3 ' )s 3 + (C 2 + C 2 ')s 2 + (C, + C/)s + Co] (5.64] 


and the zeros for velocity fluctuations are given by 

— Ra{M\s + 1 ) [. LbCbS 2 + ( R{>R c Cb + Li)s Rb + i? c ] . ( 5.65 ) 

It is significant that the terms (R C L W — R (l Lb) and ( R c Cb — R a C w ) appear in the ex- 
pressions for the transfer function coefficients given by equations (5.56) and (5.63) since 
the influence of the reactive bridge elements can be minimized if these terms approach 
zero. For both of these terms to be zero requires that Lb/L w = R c /R a = (cross-bridge 
ratio) and that Cb/C w — R a /R c = l/(cross-bridge ratio). This means that the only way 
to nullify both of these terms is to make the cross-bridge ratio unity. This may be the 
reason why a symmetrical bridge is offered as an option for high frequency operation with 
one commercially available system. The anemometer described in this memorandum can 
be configured with a symmetrical bridge as an option. 


5.5 Discussion 

The expansions given above illustrate that the derivation of transfer functions for even 
high order systems would be extremely tedious, particularly for systems with large num- 
bers of reactive bridge elements. One way to avoid the tedious algebra is to derive the 
transfer functions numerically. In essence, the numerical expressions for the system trans- 
fer functions can be built up in stages from the bridge impedances and amplifier transfer 
functions using routines for adding and multiplying the polynomials in The roots of the 
transfer function polynomials can then be solved numerically. One disadvantage of the nu- 
merical approach is that the relationships between the parameters forming the coefficients 
of the transfer function polynomials are hidden. However the analytical expressions for 
the polynomial coefficients for higher-order systems approach such a complexity that little 
can be inferred by direct examination anyway. 

It has been the authors experience that the bare minimum configuration capable of 
reproducing the behavior observed in real systems is that defined in section 5.4.3 which 
includes the frequency response characteristics of the feedback amplifier and the lumped 
inductor representing the probe cable and the balance inductor. Examples of system 
behavior are given in section 6. 
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6. EXAMPLES OF SYSTEM BEHAVIOR 


6.1 Interpretation of Electronic Square- Wave Tests 

It has been proved (see equation 5.21) that the poles appearing in the transfer func- 
tions for both offset and velocity fluctuations are identical, independent of the number of 
amplifiers in the cascade, the complexity of the overall transfer function for the feedback 
amplifier and the number of lumped components used to model the bridge impedances. 
This is extremely fortunate since the frequency response of the system to velocity fluctu- 
ations can be estimated using electronic tests. 

It is emphasized here again that one of the additional zeros in the system transfer 
function for offset voltage perturbations is always a simple zero with a time constant close 
to that of the wire filament e.g. see equation (5.37). The time constants of the system poles 
are usually at least an order of magnitude smaller than this zero. Consequently the d.c. 
component of the system response to a step input of offset voltage is greatly diminished 
with respect to the higher frequency components. Therefore this simple zero is responsible 
for causing the electronic square wave response to have an appearance resembling that 
expected from a delta-function impulse of velocity. 

The most commonly used method for tuning anemometers is to adjust the square wave 
response such that it exhibits optimally damped second-order behavior and then to es- 
timate the frequency response from the ringing frequency. However it is important to 
note that the ringing frequency is always less than the system frequency response. The 
time constant of the system is given by r 0 =2ir/\cr + ju\ while the time constant for the 
ringing is r r =27r/cj. The step (square-wave) response of a system possessing optimally 
damped complex conjugate dominant poles is shown in figure 6.1(a). Optimum damping 
for complex conjugate poles given by T 2 s 2 + + 1 occurs for £ = 0.6. For the case 

of optimum damping r 0 / r r =0.8. When the system is slightly underdamped, as shown 
in figure 6.1(b), the ringing frequency provides a more accurate estimation of the system 
frequency response. In this case, Q = 0.42, and r 0 /r r =0.91. 

While the square-wave test is an invaluable aid for tuning anemometers it is easy to 
misinterpret the actual frequency response since the system can be strongly affected by 
the higher-poles in addition to the lower-order poles appearing in the transfer function. The 
characteristic frequency of a pole is equal to the scalar distance of the pole from the origin 
of the s-plane. When there is more than one pole in the system the frequency response or 
“roll-off” frequency of the system will be defined to be equal to the lowest characteristic 
frequency of all the system poles. The pole(s) possessing this characteristic frequency will 
be referred to here as the dominant pole(s). Although the system will respond to higher 
frequency inputs it is strictly unusable beyond this point since the system sensitivity varies 
with frequency. Of particular concern is whenever a simple pole possesses a characteristic 
frequency that is lower than the complex conjugate poles. Under these conditions the 
square-wave response may still exhibit a high frequency ring but this will be at a higher 
frequency than the system response which is dominated by the simple pole. 

The results of computer simulations are used below to demonstrate the danger of blindly 
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FIGURE 6.1 Response of system with complex cojugate dominant poles. r 0 and r r are the 
time constants of the system and for the ringing, (a) Optimally damped (£ = 0.6). (b) 
Underdamped (£ = 0.42). 
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following the simplistic recipe described above. Computer simulations are a convenient tool 
for exploring the behavior of hot-wire systems since the trajectories of the poles and zeros 
on the plane, the Bode diagrams and the square-wave response can all be calculated 
for each operating point as the anemometer controls are adjusted. In a real system the 
operator can only view the square- wave response. A computer simulation allows operating 
conditions to be examined that are beyond current technology and extreme operating 
conditions can be explored without the penalty of wire breakage. In addition to illustrating 
these undesirable modes of operation, instructions are given on how to avoid them by 
appropriate adjustment of the anemometer controls. 

An interactive program has been written by the author for examining the stability and 
frequency response of constant temperature hot-wire systems. All the parameters are menu 
selectable including the number and type of components forming the impedance of each arm 
of the bridge and the fluid and wire properties. The poles and zeros of the two equivalent 
amplifiers for the feedback can be specified and the transfer functions are calculated using 
equations (5.16) and (5.18). For a given set of parameters, the static operating point is 
obtained by iteration using equations (5.4) and (5.5) allowing the lumped time constant 
of the wire T w to be evaluated. The impedance of each arm of the bridge is calculated 
as the ratio of two polynomials as expressed in equations 5. 13(a). ..(d). The coefficients of 
Zi(s), Z 2 {s) and Z$(s) given by equations (5.22), (5.23) and (5.28) are calculated using 
general purpose routines for multiplying and adding polynomials. These routines are then 
used to calculate the coefficients of the polynomial expressions for the transfer function 
zeros Q e {s) and Q u ($) using equations (5.20) and (5.27) and the poles P(s) using equation 
(5.21). The roots of the transfer function polynomials (i.e. the poles and zeros) are found 
numerically using the Lin-Bairstow method of successive quadratic factorization. If the 
system input (i.e. u or e 3 ) is sinusoidal then the output will also be sinusoidal after a 
sufficient time has elapsed for the transients to have decayed. This corresponds to putting 
s — ju in the transfer function which gives a complex number for t 0 /e s and e 0 /u . The 
absolute values \e 0 /e s \ and |e 0 /u| give the Bode diagrams for amplitude. The response 
to a more general class of inputs is given in Appendix A. The square-wave response is 
calculated using equation (A1.16). 

6,2 Offset Voltage of the Feedback Amplifier 

The results of a systematic parametric study suggest that only two types of dominant 
pole s-plane trajectories are observed as the amplifier offset voltage E qi is varied. The 
type of dominant pole trajectory depends on the nature of these poles as the frequency 
response of the feedback amplifier Ja — > oc. The effect of varying E g i on the higher-order 
poles is usually very small. 

Two examples that are representative of each type of behavior have been calculated using 
the 5 th -order system model derived in section 5.4.3. The frequency response of the feedback 
amplifier is assumed to be fiat through to a simple second-order roll-off i.e. M\ — 0 and 
M 2 = M 3 = M in equation (5.51). In each example the gain I\ and frequency response /a 
of the feedback amplifier are fixed at K = 1000 and /a ~ 79.6 kHz (i.e. M — 2 x 10~ 6 sec) 
and the trajectories of the system poles on the s-plane are traced out as the offset voltage 
Egi is varied. Only quadrant 2 of the s-plane is shown since quadrant 3 is the mirror image 
of 2 being reflected about the real axis. 
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FIGURE 6.2 Calculated s-plane trajectories (2 nd -quadrant only) of dominant poles for 
5 th -order system with increasing offset voltage E ql . Inductor values L b = 8pH, L w = 1/iH, 
amplifier frequency response f A sa 79.6kHz and gain A'=1000 are constant. Flow velocity 
U=20 m/s. Point f A = cx corresponds to model of Perry and Morrison (i.e. M = 0) for 
conditions at point 1. Square- wave response for operating points 1 to 5 also shown. Line 
C=0.6 shows optimum damping for complex poles. 
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FIGURE 6.3 Bode diagrams and step-response for point 5 in figure 6.2 where E q% =100mV, 
(a) and (b) offset voltage perturbations e a , (c) and (d) for velocity fluctuations u. 

Initially, in the first example, at the point labeled 1 in figure 6.2, the system is dominated 
by a pair of complex conjugate poles labeled A (with image B). There is also a simple pole 
labeled C and another two higher-order simple poles D and E which remain well beyond 
the frequency range of interest and their trajectories are not shown. Note that if the 
frequency response of the feedback amplifier were increased to infinity while the other 
parameters were held fixed the system would still be dominated by a pair of complex 
conjugate poles but they would be more damped as illustrated by point labeled f A = oc 
in figure 6.2. Increased damping of the dominant poles with increased f A appears to 
be a general property of constant temperature hot-wire systems. This rather surprising 
property will be discussed in more detail in section 6.4. Returning to conditions at point 
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1, as E qi is increased, the complex conjugate poles A and B move towards the real axis 
where they eventually meet and split to form two simple poles. These points are known as 
breakaway points in the system theory literature. As E qt is further increased the simple 
pole C, which has been moving towards the origin, now merges with the simple pole B to 
form a new complex conjugate pair. However pole A has continued to move closer to the 
origin and therefore it dominates the system frequency response. This behavior is typical 
of systems in which the dominant poles remain complex as /a -> oo. 

The square-w r ave response corresponding to the points labeled 1 to 5 are also shown in 
the figure. The square- wave response can be very misleading under these conditions. For 
example, for conditions corresponding to point 5 in figure 6.2, where E q , = 100roV\ the 
ringing frequency is around 320 kHz, but the simple pole A limits the system frequency 
response to about 3.6 kHz. The Bode diagrams and the step-response for both offset 
voltage perturbations and velocity fluctuations corresponding to this situation are shown 
in figure 6.3 and they clearly demonstrate why this type of system behavior is undesirable. 

An example of the second type of trajectory is shown in figure 6.4. For conditions 
corresponding to point 1 the dominant poles are also complex conjugates. However in the 
limit /a — * oo this system possesses ALL simple poles and different s- plane trajectories 
are observed as the offset voltage E q i is varied. With increasing E q i poles A and B 
remain complex conjugate while simple pole C moves towards to the origin and eventually 
dominates the system. This behavior is typical of systems which possess only simple poles 
when /a — ► oo. Despite the different trajectories of the dominant poles the end result is 
much the same as in the first example and the system possesses similar undesirable Bode 
diagrams and step-response characteristics. 

It turns out that for the 3 rd -order system considered by Perry and Morrison (1971) (i.e. 
/A = oo, see section 5.4.2), the type of behavior observed in the above examples can only be 
produced when the balance inductor Lf, is set to be very close to the value required for a.c. 
bridge balance i.e. Lbo = (72 c / R a )L w . For values of Lb slightly less than this, the system 
is dominated by the usual complex conjugate pole pair. For values of Lb slightly greater 
than Lbo the poles rapidly cross over to the RHS of the s-plane and the system becomes 
unstable. However in a real system the phenomena described above can be produced over 
a wide range of balance inductor settings by increasing the offset voltage E q i. The only 
way to simulate these observations is to include the frequency response characteristics of 
the feedback amplifier. 

6.3 Balance Inductor 

As mentioned above, the analysis of Perry and Morrison predicts instability when the 
balance inductor Lb is in excess of the value required for a.c. bridge balance Lbo ■ However 
systems with finite frequency response amplifiers are capable of maintaining system sta- 
bility when the balance inductor exceeds the a.c. balance value Lbo- For a fixed amplifier 
frequency response /a it is possible to obtain significant improvements in the system fre- 
quency response f a by suitably adjusting the offset voltage as Lb is increased to values in 
excess of Lbo- In fact having control of only E q , and Lb requires that Lb > Lbo to obtain 
the optimum system response. However if Lb is too large it may be impossible to obtain 
a satisfactory response by adjusting the offset voltage. 
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FIGURE 6.4. Same as figure 6.2 but with Lb = 0.5/^H and L w — 0.1/iH. Despite the 
different 5-plane trajectories the end result is much the same as in figure 6.2 for large E qi . 

Lb/Lbo =0.5 in the example shown in figure 6.4 and this configuration will serve as a 
baseline from which to illustrate these effects. The 5-plane trajectories for this configura- 
tion are reproduced in figure 6.5 along with those of five other configurations in which the 
values of Lb/Lbo are held fixed at 1, 2, 3, 4, and 5 respectively. For small values of the 
offset voltage (e.g. E qi < 2 mV) nearly all the configurations are unstable since the poles 
A (and B) are located on the right-half-plane. As E g , is increased these poles cross to the 
left- half- plane and the simple pole C of each system moves along the real axis towards the 
origin. For a given offset voltage, the initial position of pole C is further from the origin 
for smaller values of Lb- However this effect is pronounced only when E qi is small. As 
Eqi is increased the location of the simple poles C become more and more independent of 
Lb . For each Lb, the best response will be defined to be when the characteristic frequency 
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FIGURE 6.5. Same as figure 6.4 where L w — 0.1/iH and Lbo = (R c / R a )L w = 1/iH. 
Trajectories for Lb /Lb o= 0.5, 1, 2, 3, 4 and 5 also shown. Dashed line shows locus of 
points where poles A, (B) and C have equal characteristic frequencies. Having E g i and Lb 
as the only tunable parameters leads to an optimum response when Lb ~ l.STto* 

(i.e. | a + ju>\/2n) of poles A and B and simple pole C are equal. For Lb < Lbo the best 
response occurs when the poles A (and B) are overdamped. When Lb > Lbo the poles 
A (and B) never reach the optimum damping condition so that the best response occurs 
when they axe underdamped. The locus of the operating point with best response that can 
be obtained by adjusting both the offset voltage E g i and the balance inductor Lb is shown 
as a dashed curve. Optimum damping ((=0.6) occurs when Lb ~ 1.5L&0* The square- wave 
response corresponding to various operating points axe shown in figure 6.6 and these may 
help users interpret anemometer behavior. In paxticular note how there is no adjustment 
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of E q , that will increase the damping of system to a satisfactory level when Lb is very 
large. An operator observing this type of square- wave response is advised to reduce the 
size of the balance inductor. 


6.4 Frequency Response and Gain of the Feedback Amplifier 

Intuitively it might be expected that the roll-off frequency of the feedback amplifier 
would play a dominant role in determining the overall system frequency response. Anal- 
ysis predicts that the frequency response of hot-wire systems /„ — > oo in the limits of 
/a — ► oo, K — *• oo, E qi — > 0 and Lb — ► Lbo- However, attempts to approach these limits in 
practice axe invariably frustrated by the appearance of instabilities. For example, Smits 
and Perry (1980) observed that hot-wire systems are prone to instabilities as Lb — ► Lbo 
since there is an extreme sensitivity to very small variations in Lb- However this obser- 
vation was made for 3-order systems where / a = oo. Systems with finite frequency 
response amplifiers have an optimum frequency response when Lb > Lbo and stability can 
be maintained even when Lb is excessively large. Nevertheless one might suspect that the 
form of dominant pole instability described by Smits and Perry would eventually occur 
as /a is increased. However the results of a systematic parametric study suggests that 
other higher-order instabilities are more likely to arise beforehand. While there are still 
significant improvements to be found by increasing the frequency response and gain of the 
amplifier, its role is often more critical in determining the system stability by influencing 
the higher-order poles. Two examples are given in figures 6.7 and 6.8. 

Calculated s-plane trajectories of a typical 5 th -order system are shown in figure 6.7. 
Two cases are considered i.e. increasing amplifier frequency response and increasing gain. 
The behavior of the dominant poles labeled A and B and the higher-order poles labeled 
C, D and E are shown separately in figures 6.7 (a) and (b) because of the wide range of 
values observed. In the first case, the amplifier gain A'=1000, the offset voltage £,,=12.5 
mV and the inductors Lb = 40pH and L w = 5pH are held constant. Note that the 
damping of the dominant poles A and B increases with increasing amplifier frequency 
response. As a matter of interest the system frequency response (f 0 ssl7.8 kHz) is higher 
than the amplifier frequency response (/ a ~15.9 kHz i.e. M = 1.0 x 10 -5 s) for conditions 
corresponding to point 1A labeled in figure 6.7. The system response increases rapidly 
with increasing amplifier frequency response but only to f a «27 kHz for Ja sj 106 kHz 
(i.e. M = 1.5 x 10~ € sec). Further increases in the amplifier frequency response only have 
a small effect on the dominant poles e.g. f 0 — >31 kHz for /a — ► oo (i.e. M — > 0). As 
the amplifier frequency response is increased two of the higher-order simple poles C and 
D move towards each other and merge to become complex conjugates for / a ~ 106 kHz. 
This point is labeled 2A in figure 6.7(b). The third simple pole labeled E moves further 
away from the origin and exerts negligible influence on the system. With further increases 
in the amplifier frequency response the two higher-order conjugate poles C and D move 
towards the imaginary axis and ultimately cross over to the RHS of the s-plane resulting 
in instability. 

Similar behavior is observed if the amplifier gain K is increased while the other pa- 
rameters are held constant as shown in figure 6.7(b). Although the dominant poles move 
along slightly different trajectories the endpoints K — oo and f \ = oo are very closely the 
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FIGURE 6.7. Calculated 5-plane trajectories of system poles for increasing / 4 (with 
constant A'=1000) and for increasing K (with constant % 79.6 kHz). E qt = 12.5mV, 
L w — 5^H and Lb — 40pH. (a) Dominant poles (b) Higher order poles. 
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FIGURE 6.8. Same system as figure 6.6 but E q% increased to 20mV. Calculated s-plane 
trajectories of system poles for increasing f & (with constant A”=1000) and for increasing 
K (with constant ~ 79.6 kHz), (a) Dominant poles (b) Higher order poles. 

same. The s-plane trajectories of the higher-order system poles C and D are very similar 
as A — ► oc and /^ — > oc. 

In the previous example the dominant poles approach the complex conjugate 3 rd -order 
system values asymptotically as K — + oo and as / .4 — ► 00. In the second example shown 
in figure 6.8(a), a different type of trajectory is observed i.e. the dominant poles A and 
B move towards the real axis where they meet and split to form two simple poles. These 
poles also approach the 3 rd -order system values asymptotically as AT — ► 00 and as f.4 — + oc. 
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This type of behavior is also consistent with the observation mentioned previously that 
relates the increased damping of the dominant poles with increased frequency response 
and gain of the feedback amplifier. Note that increasing /a or A beyond the point where 
the dominant complex poles are transformed into two simple poles actually reduces the 
overall system frequency response since pole A moves closer to the origin. The effects of 
increasing A and / a on the higher-order poles are very similar to those in the previous 
example, despite the different trajectories for the dominant poles, i.e. the poles C and D 
merge to become complex conjugate and move towards the imaginary axis. The higher- 
order poles ultimately cross over to the RHS of the 5-plane resulting in system instability 
as shown in figure 6.8(b). 

These examples demonstrate several important properties of constant temperature hot- 
wire systems. The effects of increasing A and / a are quite similar and have the unexpected 
effect of increasing the stability of the dominant poles. This appears to be true for all sys- 
tems i.e. the author has not found one contrary example to this observation out of a large 
number of cases studied. It is the higher-order poles that are responsible for the system 
instability under these conditions. It is often difficult to determine the damping of the 
higher-order poles during a square- wave test, even when they are grossly underdamped, 
since the oscillations can still decay rapidly when compared to the response of the dom- 
inant poles. However a small change in either K or /a from this point could result in 
instability. This type of higher-order instability can occur suddenly and without warning 
to the anemometer operator who can only observe the square- wave response. 

6.5 Instabilities with Subminiature Wires 

Some workers (e.g. Miller, Shah and Antonia 1987) have reported frustration with 
instabilities when using subminiature hot-wires. One way of comparing system behavior 
with different diameter wires is to non-dimensionalize the poles and zeros with the lumped 
time constant of the wire T w . Alternatively, the coefficients of s n in the system transfer 
functions can be made non-dimensional by division with R w 2 T w n . The relative positions 
of the non-dimensional poles and zeros will be the same for systems in which R a /R w , 
Rb/Rw , Rc/Rw , R/Rw, L b /(R W T W ), L W /{R W T W ) and M/T w are equivalent. 

For a given l/d ratio, R g cx 1/d, so that R g for a typical subminiature wire (e.g. d=0.5 
fim) is about an order of magnitude greater than a more conventional wire (e.g. d=5 
/on). T w is reduced for smaller diameter wires, but by a smaller factor than the length- 
scale ratio since the Reynolds number is also reduced. Consequently the product R W T W 
becomes larger as the wire is made smaller. This leads to a beneficial reduction in the size 
of the non-dimensional bridge inductance Lb/{R W T W ) and L W /(R W T W ). Values of R a /R w , 
Rb/R w , Rc/Rw and R/R w can be altered by adjusting the bridge resistors and the offset 
voltage. Finally, the frequency response of the feedback amplifier must be increased to 
obtain the same value of M/T w . 

If /a is fixed then M/T w can be considerably larger with subminiature wires. This is 
especially significant since it has been demonstrated that the dominant poles become less 
damped as f \ is reduced. If /a is too small then the dominant poles may be unstable. 
Figure 6.9 shows the effect of reducing the wire diameter while maintaining the same 
length-to-diameter ratio (1:200) and the same resistance ratio R « 2 (by increasing R b ). 
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FIGURE 6.9. Calculated 5-plane trajectories of the dominant poles for a range of wire 
diameters with the same length to diameter ratio. Rb has been adjusted to closely give the 
same resistance ratio i.e. R w 2. E qt = 5mV, 1,6=0. 5 /z H and L w —0.1 /zH. Higher values 
of K and Ja are required for stability as the wire diameter is reduced, (a) Increasing / a 
(with constant gain A'=1000). 

At the lowest values of /a and K the dominant poles are unstable for small wire diameters. 
As either / a or K is increased these poles move to the LHS of the s-plane and the systems 
become stable. The effect on the higher-order poles is similar to the examples shown in 
figures 6.7 and 6.8. 

The wire current required for a given resistance ratio is smaller for subminiature wires 
so that the static output voltage at the top of the bridge is considerably less than that 
obtained with conventional wires. Therefore increasing R a and R c to obtain the same 
values of R a /R w and R c /R w , as suggested above, would also help restore the size of the 
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FIGURE 6.9. Continued, (b) Increasing K (with constant /a ~ 79.6 kHz). 

output signal. However this is not recommended since increasing R a and R c tends to reduce 
the damping of the dominant poles as shown in figure 6.10. The largest values of R a /R w 
and R c /R w correspond to the values used for the d=5.0/L/m wires in earlier sections. An 
even higher frequency response amplifier would be required to obtain a satisfactory system 
response under these conditions. 

The author suggests that insufficient frequency response of the feedback amplifier is the 
most likely cause of system instability when using subminiature wires. 

6.6 Effects of Bridge Capacitance 

The previous examples used a 5 th -order system in which the frequency response of the 
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FIGURE 6.10. Increasing upper bridge resistors R a and R c causes dominant poles to 
become less damped. Subminiature platinum filament of length 0.1mm and diameter 
d=0.5/wn ( R g ss 8017). Largest values of R a /R w and R c /R w axe approximately the same 
as for the d=5//m wires in figures 6.2 to 6 . 8 . Lj> = 8 /uH, L w = 1 //H, Rb = l.QkQ. and the 
amplifier frequency response /a ~ 1 . 6 MHz, gain K =1000 and offset voltage E qx = 5mV 
are constant. Air velocity is 20 m/s. 

feedback amplifier is assumed to be flat through to a simple second-order roll-off i.e. M\ = 0 
and M 2 = M 3 = M in equation (44). In this section an example of the effect of bridge 
capacitance will be presented. The transfer functions for a system with capacitance in the 
active and balance arms of the bridge wets expanded in section 5.4.4. 

Figure 6 . 11 (a) shows the step response and Bode diagram for offset voltage perturbations 
corresponding to the 5 th -order system defined in figure 6.7 (i.e. Cb = C w = 0), where 
the gain and the frequency response of the feedback amplifier are A ’ = 1000 and f a =477 
kHz (i.e. M = 3.33 x 10 ' sec). The spike appearing in the Bode diagram at around 
5 MHz is the result of the higher-order complex conjugate poles which are close to the 
imaginary axis. Although these poles are very much underdamped, the oscillations caused 
by exciting the system with a square- wave decay rapidly when compared to the response of 
the dominant poles. When a small capacitor ( C w = 5 pF) is placed in the active arm of the 
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bridge these poles move even closer to the imaginary axis while the effect on the dominant 
poles is barely detectable. An additional simple pole is introduced since the system transfer 
function is now 6 th -order but it is located well away from the region of interest and exerts 
very little influence on the system. As a result of the movement of the higher-order poles, 
the spike in the Bode diagram is much sharper as shown in figure 6.11(b). The envelope 
for the exponential decay of the oscillations is now much slower than the dominant poles 
so that high frequency oscillations are clearly visible in the step-response as shown in 
figure 6.11(b). If a small capacitor is added to the balance arm of the bridge (C& =10pF) 
the system transfer function becomes 7 th -order. The higher-order simple pole originally 
present in the 5 th -order system forms a complex conjugate with the additional pole but 
this pair is close to optimum damping and they are located far away from the region of 
interest. The simple pole resulting from the addition of Cu,=5pF is unaffected by the 
addition of Cj and once again the effect on the dominant poles is negligible. However the 
complex conjugate poles which were moved closer to the imaginary axis are now moved 
away towards the optimum damping condition resulting in almost complete removal of the 
high frequency spike in the Bode diagram as shown in figure 6.11(c). This example serves 
to demonstrate that even small quantities of bridge capacitance can have a dramatic effect 
on the system stability but very little influence on the dominant poles which determine 
the overall frequency response of the system. 
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APPENDIX A. LAPLACE TRANSFORM METHODS 


Al.l Introduction 

The hot-wire anemometer is a difficult instrument to understand properly. Subtle in- 
sights into the system behavior are often required to detect an erroneous mode of operation. 
A major source of difficulty is that the electronic square-wave test is an indirect method 
for estimating the frequency response. While the poles of the transfer function for velocity 
and offset voltage fluctuations are identical, the zeros are quite different and this is often 
the source of considerable confusion. 

Once the transfer function is known, the behavior of complicated systems of differential 
equations can be analyzed in the so called frequency domain by using simple algebra. By 
taking inverse transforms, the complete time domain solutions can be obtained by applying 
the simple algorithms that are described below. These techniques have been adapted into 
a computer simulation which is used to examine how models of the Constant Temperature 
Hot-wire Anemometer respond to various forcing conditions. The advantage of a computer 
simulation is that the trajectories of the poles and zeros on the 5-plane, the Bode diagrams 
and the square-wave response, can all be calculated for both velocity and offset voltage 
fluctuations as the anemometer controls are adjusted. With a real system the operator 
can only observe the electronic square-wave response which is open to misinterpretation. 


A1.2 Definitions 

Given a function of time f(t ), then the Laplace Transform L{/(t)} is defined as 


M /(<)> = / me-'dt = f{s) (.4.1) 

Jo 

The Laplace variable s is in general a complex quantity i.e. s = a + juj. The Laplace 
transform of a derivative has the extremely useful property that 


L 



= s n f(s) 



+ 5 


n-2 4f(<) 

dt 


+ -•• + 


t= o 


dt n ~ l 



(A.2) 


For systems in which all initial conditions are zero, the transfer function can be derived 
by expressing the differential equation in terms of the input and output variables and 
simply replacing the differential operator d n f(t)/dt n with s n . The equation can then be 
manipulated by treating s as an ordinary algebraic quantity. The following differential 
equation is used to provide some definitions. 


A n 


d n e 0 

dt n 


+ A n -i 
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dt n ~ l 
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+ B m - 1 
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dt 171 - 1 


-(-••• + e, 


(.4.3) 


The input e* and its derivatives on the r.h.s side of (A. 3) are known as the forcing function. 
Assuming that all initial conditions are zero, then applying equation (A.2) gives 


e ° t \ r- + B m -\S m 1 + ' • ’ + 1 

e; v A n s n + An-is”- 1 + ••• + ! 


(AA) 
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When dealing with electrical circuits the concept of operational impedance allows the 
transfer function to he constructed without having to write down the differential equations. 
The operational impedance of a capacitor is 1 jCs. The operational impedance of an 
inductor Ls. The transfer function for a circuit (i/i) can then be evaluated by using 
Kirchoff's Laws. 

Some alternative forms of the transfer function are given below: 


o 4. (jW-i i . . 

^ B m 3 + 

• + TT 

. 4 „ + + " 

B m (s — z n )($ — z n — i ) - • 

• + t 

■(s~z l ) 

Art ( Pm)($ Pm — l)' 

(T Zn s + l)(T Zn _ l s + 1 ) • • 

■■ ($ - pi ) 

' • (T.j s + 1 ) 

(Tp m S + \){Tp m _ 1 5 + 1 ) • 



(.4.5) 


(.4.6) 


(.4.7) 


The quantity B m /A n is often called the level factor. The quantity zj is called a zero of 
the system since the transfer function goes to zero when s = Zj. The quantity p } is called 
a pole of the system since the transfer function becomes infinite when s = p v T Pj is the 
time constant of pole j and TV is the time constant of zero j. The poles and zeros may 
be real or complex quantities. If they are complex they must occur in conjugate pairs for 
the transfer function polynomial coefficients to be real numbers. If a pole or zero is real 
then it is referred to as a simple pole or a simple zero. 


A1.3 System response to forcing 

If the input e{ is sinusoidal then the output e a will also be sinusoidal after a sufficient 
time for the transients to decay. In this case the system response can be calculated by 
putting s = ju; where ^ is the forcing frequency in rad/s. This gives a complex number 
for (e 0 /e,-)(w). A plot of the amplitude versus frequency is known as the Bode diagram. 
The response to a more general class of inputs was considered by Perry (1982) using 

e t =E t H(t)e Xt (AS) 

where E, and A are in general complex quantities. The Heaviside step function H(t) 
ensures that initial conditions are zero. The Laplace transform of the input is given by 


Therefore the output of the system described by (A.4) is given by 

. i x KE t B(s) 

Co{s) = -T-j 7 — r; TV (A9) 

pn ) ’’’(*§ Pi)(«s A) 

where B(s) = B m s m + B m - i s rn ~ 1 + • • • + 1. Once the denominator has been factorized 
the equation can be expressed as partial fractions, 
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The coefficients C„ ■ ■ • C\ and C\ can be found by equating coefficients of ( A. 9) and (A. 10). 
By allowing s to approach p„ ,p „-] , • • ■ p\ and p\ in turn then the appropriate term in 
(A. 10) will dominate. Hence 


= < An > 

C A = Km{EiF{s)} (.4.12) 

where F(s) is the full transfer function as given in (A. 4) to (A. 7). The partial fractions 
in equation (A. 10) allow the inverse transform to be obtained rather simply so that the 
output can be expressed as a function of time 


e 0 (f) 

Ei 



F(s)( s - Pj n t 

(s-X) } 


+ F(\)e Xt 


(.4.13) 


It can be seen that the solution consists of two parts. The first part is the complimentary 
function and the time constants depend on the poles. The second part is the particular 
integral. The response rate depends on A which is a property of the input. Equation (A. 13) 
is a slightly generalized version of the Heaviside expansion theorem and is given by Perry 
(1982). 

The full transient response of the system to the sinusoid 


e,(t) = EiH{t)e juJ 1 


(A. 14) 


can be found from (A. 13) by putting A = ju;. The square- wave response he. the response 
to the step input 

€i(t) = EiH(t) (A. 15) 

can be found by putting A = 0 which leads to 

^ £ lim (£^( a - p> )}cW f + F(0). (A.16) 

E, l s J 

The response to a delta function e, = EsS(t) can be found by writing this as 

e,- = EsH(t) lim /——l (A.17) 

9— 0 l 9 J 


which gives 


= it, lim { F ( 5 X 5 ~Pj)) ePit + /g)E(-l/q)e t/q 

. 7 = 1 


(.4.18) 
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